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Abstract  ( cont . ) 

The  NORSAR  Deccection  Processing  System  has  been  operated  throughout 
the  reporting  period  with  an  average  u[ititne  of  98.2  per  cent.  A  total 
of  1958  seismic  events  has  tieen  reported  in  the  NORSAR  monthly  seismic 
bulletin.  The  performance  of  the  continuous  alarm  system  and  the 
automatic  bulletin  transfer  by  telex  to  AFTAC  has  been  satisfactory. 
Processing  of  requests  for  fvill  NORSAR/NORESS  d.ita  on  magnetic  tapes 
has  progressed  according  to  estal)li.shed  schedules. 

The  satejlite  link  for  transmitting  NORESS  data  in  real  time  to  the 
U.S.  has  had  an  average  uptime  of  99.3  per  cent.  On-line  NORESS 
detection  processing  and  data  recording  at  the  NORSAR  Data  Center 
(NDPC)  has  been  conducted  throughout  the  period,  with  an  average 
uptime  of  97.8  per  cent. 

The  ARCESS  array  started  operation  in  mid-October  1987,  and  the  data 
v;ere  initially  recordf?d  and  processed  at  NDPC  using  a  Sun  2-based 
computer  system.  In  May/June  1988,  this  system  experienced  signif  i  c.ant 
hardware  problems,  and  the  time  schedule  for  the  planned  edaange  to  .i 
Sun  3  system  wa.s  therefore  .iccelerated.  The  changeover  was  successful  I  v 
completed  in  early  July.  Average  recording  time  for  ARCESS  was  77.2% 
for  the  total  report  i;:g  period,  a.vl  per  cent  when  disregarding  tht 
May 'June  period 

r'iild  maintenance  activity  lias  iiicludid  regular  pre'.'entive  maintenatice 
at  ;:11  array  si  te  r-  and  occasional  correct  ive  actions  when  required.  In 
;',;r '  icular ,  much  work  has  been  conducted  at  the  ARCESS  array  site, 
including  instillation  of  a  new  Clob.-i]  Positioning  System  synchronlzt  d 
clock  in  coope.r.it  ion  with  Saiv.lia  personnel  and  removal  of  the  KS-3600‘': 
borehole  se  i  .smoineter  ;or  repair.  The  .NDRSAR  and  NORESS  field  sysCom.s 
performed  entlrcdy  s.i  t  i  s  Lac  t  or  i  1  v  ibronghout  the  reporting  period. 

A  considerable  etfort  has  been  oxponded  in  upgrading  the  oti-line  and 
off -line  de tec t ion/evi nt  processing,  software  which  is  being  developed 
at  NORSAR  for  general  array  applications.  The  program  systems  have  ht:er, 
tested  on  data  from  NORS.AR,  NORESS,  ARCESS,  FINESA  and  Crafenberg,  ami 
the  implementation  will  lie  toordinated  with  the  Intelligent  Array 
S  vs  tern  deve  lopirei:  t  .s  . 

.A  study  of  Lg  spi'Ctra  of  NORSAR  -  rceorded  exp  Ic.s  i  on.s  from  the  Shagan 
River  ter,',  arta  near  Semi  p.',  1  .i  t  i  n.'-'r  ,  USSR,  has  sliown  :h;it  the  main 
energy  in  the  Lg  w<'ivetrain  is  conlined  to  the  fr('(|nency  range  0 , 6  • 

1.0  Hz.  There  is  some  rvidence  of  source  size  scaling  effects  (i .c  . . 

dominant  Lg,  f  I't'iiuencv  for  l,irp,er  events),  but  the  variation  is 
.small  and  apjH'.ir.s  to  'oe  of  litt  Ic  ■;  i  gu  i  f  i  cance  in  RMS  Lg  magnitude 
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estimation.  The  Lg  spectra  show  no  significant  differences  for  events 
from  the  two  portions  of  the  test  site  (NE  Shagan  and  SW  Shagan) .  On 
the  other  hand,  the  spectral  difference  between  early  P  coda  and  Lg  is 
larger  for  SW  events  than  for  NE  events,  and  this  holds  true  in  the 
entire  frequency  band  0.6 -3.0  Hz. 

An  analysis  has  been  made  of  statistics  on  P,  PcP,  PKP  and  PKKP  travel 
time  residuals  with  respect  to  the  isotropic  PREM  earth  model  using  ISC 
bulletin  data  for  the  years  1964-1984  as  well  as  other  data  sources. 

The  scatter  in  the  residuals  is  significantly  greater  for  shallow 
events  than  for  deep  events.  For  P  phases  at  distances  less  than  65° 
this  increased  scatter  reflects  increases  in  number  of  "early"  as  well 
as  "late"  arrivals,  whereas  for  more  distant  P  phases  and  other  phases 
the  increase  in  scatter  is  dominated  by  a  larger  number  of  late 
arrivals.  The  statistics  further  show  that  after  applying  the  isotropic 
PREM  model  there  is  still  a  significant  mean  residual  left  in  the  data, 
and  this  is  particularly  pronounced  for  core  phases. 

The  coupling  mode  technique  for  modelling  surface  wave  propagation  in 
2-D  structures  presented  in  previous  Semiannual  Technical  Summaries  has 
been  applied  to  a  model  of  the  North  Sea  Graben.  The  purpose  has  been 
to  examine  how  a  large-scale  and  very  strong  lateral  variation  of  the 
crustal  structure  affects  the  propagation  of  the  short-period  surface 
wave t rains . \ Our  continued  work  has  confirmed  the  conclusions  given  in 
earlier  repotrts.  Thus,  our  numerical  modelling  of  Lg  wave  propagatioti 
in  a  simplified  model  of  the  North  Sea  Central  Graben  does  not  predict 
the  severe  attenuation  of  the  wavetrain  actually  observed  in  this 
region.  On  the  contrary,  the  Lg  wavetrain  appears  very  robust  when 
crossing  a  zone  where  its  waveguide  is  strongly  deformed.  Since  the 
large-scale  geometry  of  the  Graben  fails  to  explain  the  observed  data, 
it  is  suggested  that  future  work  explore  alternative  explanations  for 
the  observed  attenuation,  such  as  scattering  by  2D  or  3D  basaltic 
intrusions  in  the  lower  crust,  extensive  faulting  associated  with 
intra-fault  weak  material,  or  more  rheological  aspects  of  this  problem. 

A  detailed  analysis  has  been  made  of  the  recent  (August  8,  1988) 
earthquake  offshore  Norway,  which  was  recorded  at  NORSAR,  NORESS, 

ARCESS  as  well  as  ^  large  number  of  stations  at  regional  and  tele- 
seismic  distances  .‘^This  earthquake  is  the  largest  in  the  region  for  at 
least  30  years,  with  an  estimated  mj,  —  5.2.  Our  focal  mechanism 
solution  indicates  thrust  faulting  along  a  NNE-SSW  striking  fault 
plane,  in  response  to  E-W  compressional  stress.  The  seismic  moment  has 
been  estimated  at  10^^  Nm,  with  indications  of  scaling  consistent  with 
an  w-square  source  model.  A  major  source  of  uncertainty  in  this 
analysis  is  tied  to  the  lack  of  accurate  knowledge  of  the  anelastic 
attenuation. 
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In  cooperation  with  the  University  of  Helsinki,  an  analysis  has  been 
conducted  of  events  recorded  by  the  three  regional  arrays  in  Fenno- 
scandia  (NORESS,  ARCESS,  FINESA) .  The  latter  array  has  been  recon¬ 
figured  to  closely  resemble  the  subgeometry  of  NORESS  and  ARCESS 
obtained  when  excluding  the  D-ring.  As  a  first  part  of  this  study,  the 
detection  capabilities  of  the  upgraded  FINESA  configuration,  using  the 
RONAPP  algorithm,  was  investigated.  The  results  were  excellent,  in  that 
98  out  of  103  reference  events  (listed  in  the  Helsinki  bulletin)  had  at 
least  one  detected  P  or  S  phase  at  FINESA.  Many  additional  events  were 
also  detected.  The  second  part  of  the  study  has  addressed  the  location 
capabilities  of  the  three-array  network,  again  using  the  Helsinki 
bulletin  as  a  reference  for  comparison.  On  the  average,  joint  three- 
array  locations  deviate  from  the  reference  location  by  only  16  km, 
whereas  comparable  deviations  using  one  or  two  arrays  for  location 
purposes  were  68  and  26  km,  respectively.  Since  the  arrival  times  used 
were  those  determined  automatically  by  on-line  processing,  there  are 
clear  possibilities  for  further  improvements  in  multi-array  location 
accuracy.  This  could  be  achieved  both  through  more  accurate  readings, 
introduction  of  more  detailed  regionalized  travel-time  tables  and 
application  of  master-event  location  techniques. 

.As  a  continuation  of  earlier  studies  on  Lg  magnitudes  and  yield 
estimation  of  Semipalatinsk  explosions,  we  have  now  completed  the 
analysis  of  all  available  Gr^fenberg  Lg  recordings  of  large  Shagan 
River  explosions.  Because  of  the  varying  number  of  channels  available 
and  the  relatively  large  systematic  variations  in  Lg  amplitudes  across 
the  array,  we  have  used  station  correction  terms  for  individual 
channels  when  computing  average  magnitudes.  The  resulting  values  show 
excellent  correspondence  with  NORSAR  Lg,  with  the  standard  deviation  of 
magnitude  differences  between  the  two  arrays  being  about  0.03  units  for 
well - recorded  events  (more  than  b  of  the  13  GRF  channels  available).  Ue 
have  also  investigated  the  effect:  of  applying  individual  instrument 
correction  terms  to  NORSAR  Lg  data  and  have  found  the  effects  to  be 
modest.  A  study  comparing  joint  NORSAR/Gra fenberg  Lg  magnitudes  to 
maximum- likelihood  mj^  based  on  ISC  data  has  confirmed  the  pattern 
previously  observed,  i.e.,  a  low  P-Lg  bias  for  NE  Shagan  in  comparison 
to  SW  Shagan.  The  JVE  explosion  of  lA  September  1988  had  a  P-Lg  bias 
close  to  the  average  for  the  SW  Shagan  event  set.  This  explosion  was 
estimated  at  m(Lg)  =  5.97  at  both  NORSAR  and  Grafenberg. 
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SUMMARY 


I  , 

This  Final  Technical  Summary  describes  the  operation,  maintenance  and 
research  activities  at  the  Norwegian  Seismic  Array  (NORSAR) ,  Norwegian 
Regional  Seismic  Array  (NORESS)  and  the  Arctic  Regional  Seismic  Array 
(ARCESS)  for  the  period  1  April  -  30  September  1988. 

The  NORSAR  Detection  Processing  System  has  been  operated  throughout  the 
reporting  period  with  an  average  uptime  of  98.2  per  cent.  A  total  of 
1968  seismic  events  have  been  reported  in  the  NORSAR  monthly  seismic 
bulletin.  The  performance  of  the  continuous  alarm  system  and  the 
automatic  bulletin  transfer  by  telex  to  AFTAC  has  been  satisfactory. 
Processing  of  requests  for  full  NORSAR/NORESS  data  on  magnetic  tapes 
has  progressed  according  to  established  schedules. 

The  satellite  link  for  transmitting  NORESS  data  in  real  time  to  the 
U.S.  has  had  an  average  uptime  of  99.3  per  cent.  On-line  NORESS 
detection  processing  and  data  recording  at  the  NORSAR  Data  Center 
(NDPC)  has  been  conducted  throughout  the  period,  with  an  average  uptime 
of  97.8  per  cent. 

The  ARCESS  array  started  operation  in  mid-October  1987,  and  the  dat,. 
w“re  iniria'’’''  recorded  •’nd  processed  at  NDPC  using  a  Sun  2 -based 
computer  system.  In  May/June  1988,  this  system  experienced  significant 
hardware  problems,  and  the  time  schedule  for  the  planned  change  to  a 
Sun  3  system  was  therefore  accelerated.  The  changeover  was  successfi.il  1  v 
completed  in  early  July.  Average  recording  time  for  ARCFSS  was  77.2? 
for  the  total  reporting  period,  and  91  per  cent  when  disregarding 
May/June  period. 

Field  maintenance  activity  has  included  regular  preventive  maintenance 
at  all  array  sites  and  occasional  corrective  actions  when  required.  In 
particular,  much  work  has  been  conducted  at  the  ARCESS  array  site, 
including  installation  of  a  new  Global  Positioning  System  (GPS) 
synchronized  clock  in  cooperation  with  Sandia  engineers  and  removal  of 
the  KS-36000  borehole  seismometer  for  repair.  The  NORSAR  and  NORESS 


field  syst..-nis’  performed  entirely  .sat  isfai  t  or  i  1  y  rlirou,yhout  the 
reporting  period. 

A  considerable  effort  has  been  expended  in  upgrading  the  on-line  and 
off-line  detection/event  processing  softw.are  vhich  is  being  developed 
at  N'ORSAR  for  general  array  applications.  The  program  systems  have  1).  , 
tested  or.  data  from  NORSAR,  NORKSS,  ARCESS,  FINESA  and  Grafenberg,  ai  d 
the  implement.stion  will  he  coordinated  with  the  Intelligent  Array 
System  developments. 

The  research  activity  is  summarized  in  Section  VII.  Section  VII. 1 
give.s  results  from  analysis  of  P  coda  and  Lg  spectra  of  Shagan  River 
explosions  recorded  at  NORSAR  and  NORESS .  Statistics  of  ISC  travel  time 
residuals  in  com.par ison  to  the  PRE.M  model  are  presented  itj  Section 
vri.2.  Section  VII. 3  reports  on  modelling  of  Lg-wave  prop.agation  across 
t.iic  Central  Graben  of  the  North  Sea,  while  observed  ground  motions  arid 
inferred  source  parnmoter.s  for  the  August  8,  1988,  More  Basin  earth 
■juake  are  discu.ssed  in  Section  VI  I. A.  An  .inalvsis  of  rep.ional  seismic 
events  using  the  NORESS/ARCESS/FINESA  .irrav.s  is  presented  in  Section 
VI  [.'j.  In  Section  VII. 6  results  of  a  comp.u'ative  analysis  of  NORSAR  and 
Gr,.i  tenberg  i.g  magniriide.s  tor  Sliay.-m  River  explosions  .na-  t  ''.'en. 


NORSAR  '-"'RRATION 


1  I  . 

1  I  .  1 _ Uftectlon  Processor  (DP)  Operation 

There  have  been  breaks  iri  the  otherwise  continuous  operation  of  the 
'TIRSAR  online  system  within  the  current  f>  month  reporting,  interval.  The 
uptime  percentage  for  tin-  period  is  '>8.?  as  compared  to  96 ,  tor  th.c- 
in'evious  period. 

Fig.  II. 1.1  and  the  accompar.y  i  ng  Table  1 1.  .  1  .  1  l:)Otli  .shovt  the  'lailv  '.)? 
downtime  tor  the  days  between  1  April  and  TO  Septeiiber  1988.  Tlte 
monthly  recording  times  and  percentages  are  given  In  T:il.'’.e  1!.1.9. 

The  breaks  can  he  grouped  .,s  follows: 


a!  Hardware  failure 

h  Stops  related  to  prograiti  work  or  error  " 

e>  Hardware  maintenance  stop-.  ! 

d)  I’owor  jtimps  and  breaks  o 

et  I'uD  error  correction  i 

li  Fommunicat  Ion  lines  9. 


The  total  downtinie  for  tin  peri-u!  wa.s  8n  Hours  and  .  :i:i:!u:e.s,  I'h*  mean¬ 
time  -  be  tween  -  fa  i  I'.ire  s  (MTKK)  w,i<;  (l..v.s  ,  as  ev>mpare<!  to  ,  i  ;  u  t  lu 

■previous  period. 


CAY 


List  or  breaks  in  op  processing  the  i^st  half-year 


I;  AY 

START 

STOP 

COKNENTS . 

9  J 

9 

56 

9 

57 

LINE  FAILURE 

94 

6 

1 

6 

2 

IINE  FAILURE 

96 

6 

1 

6 

2 

TOD  RETARED  20KP 

SI 

9 

0 

9 

2 

LINE  FAILURE 

99 

5 

57 

5 

58 

LINE  FAILURE 

101 

19 

57 

19 

58 

LINE  FAILURE 

VC4 

6 

1 

6 

2 

TOO  RETARED  lOKS 

104 

9 

10 

9 

13 

LINE  FAILURE 

1  1  1 

6 

1 

6 

2 

TOO  RETARED  IVMS 

118 

1 1 

3 

11 

10 

TOD  RETARED  lOMS 

122 

14 

40 

14 

43 

LINE  FAILURE 

126 

8 

40 

0 

41 

TOO  RETARED  20MS 

130 

1 

19 

7 

21 

LINE  FAILURE 

130 

8 

26 

8 

27 

LIME  FAILURE 

132 

6 

29 

8 

30 

LINE  FAILURE 

139 

7 

0 

7 

1 

LINE  FAILURE 

140 

5 

12 

5 

16 

LINE  FAILURE 

141 

6 

1 

6 

2 

TOD  RETARED  10K3 

142 

14 

28 

14 

31 

LINE  FAILURE 

14S 

6 

1 

6 

2 

TOO  RETARED  12MS 

149 

11 

14 

11 

23 

MOOCOMP  FAILURE 

150 

16 

26 

16 

51 

POWER  FAILURE 

165 

14 

0 

IS 

5 

MODCOMP  EAILUPC 

15  9 

S 

20 

8 

40 

DISK  FAILURE 

L-.  5 

9 

12 

9 

13 

r.iNE  FAn.uRt: 

It"’ 

1  3 

4 

24 

0 

2701  FAILURI. 

168 

0 

0 

13 

30 

270  1  FAII.URL 

;t>9 

4 

24 

19 

44 

2701  FaILUKE 

172 

6 

20 

12 

48 

2701  FAILURE 

193 

1 

56 

6 

47 

POWER  failure 

186 

6 

1 

6 

2 

TOD  RETARED  13  M 

188 

7 

29 

7 

48 

CE  KAINTENANCF  2 

195 

6 

21 

6 

22 

TOO  RETARED  2  1M:. 

19" 

18 

30 

20 

7 

POWER  FAILURE 

2C  1 

16 

30 

17 

30 

POWER  FAILURE 

201 

17 

30 

24 

0 

MOOCOMP  FAILURE 

DAV 

202 

START 

0  0 

STOP 

8 

34 

217 

9 

50 

9 

51 

218 

2 

33 

3 

31 

221 

6 

2 

6 

3 

228 

9 

29 

9 

30 

229 

13 

30 

17 

50 

2  31 

b 

33 

5 

34 

232 

7 

10 

■7 

12 

234 

10 

2 

lb 

16 

235 

6 

18 

6 

20 

235 

1  1 

25 

11 

52 

235 

12 

3 

12 

4 

236 

6 

0 

6 

3 

236 

13 

28 

13 

31 

238 

13 

20 

13 

21 

238 

15 

44 

15 

45 

239 

8 

30 

8 

31 

244 

6 

2 

6 

3 

250 

b 

1 

0 

2 

rs9 

K 

1 

6 

2 

2"4 

7 

1 

7 

2 

COfWENTS . 

MODCOMP  FAILURE 
LINE  FAILURE 
PONER  FAILURE 
TOD  RETARED  46HS 
LINE  FAILURE 
POWER  FAILURE 
LINE  FAILURE 
LINE  FAILURE 
LINE  FAILURE 
LIKE  FAILURE 
LINE  FAILURE 
LINE  FAILURE 
TOD  RETARED  I 2 NS 
LINE  FAILURE 
LINE  FAILURE 
LINE  FAILURE 
LINE  FAILURE 
TOD  RETARED  15MS 
TOO  RETARED  UKS 
TOD  RETARED  22KS 
TOD  RETARED  14MS 


T.ih  1  e  1 1 .  1  .  ]  ■  Dali)  DP  dov/nt  imi 
1988. 


t  In-  per  i  od  I  A| •  ?  1  i 


Ui  Si'p  ‘  r 


Month 


DP  Uptime  DP  Uptime  No.  of  DP  No.  of  Days  DP  MTBF''-' 
hours  %  Breaks  with  Breaks  (day.s) 


APR  88 

719.41 

99 . 9 

MAY 

741.08 

99 . 6 

JUN 

672 . 20 

93.4 

JUL 

721 .07 

98.9 

.-my: 

/37  .  4.' 

99,2 

SEP 

719. 

100.0 

10 

9 

2.8 

12 

11 

2.4 

6 

6 

4 .0 

G 

7 

4 . 3 

17 

13 

1.7 

3 

j 

7 . 5 

*  Mean- 1  ime -be  tv.’een  -  fai  lutes  total  uptime/no.  of  up  intervals. 


O’.'.line  system  performance,  1  April  - 


Table  II. 1.2. 


30  September  1988 


II.2 


Array  communications 


Table  11.2.2  reflects  the  performance  of  the  system  through  the 
reporting  period.  All  the  subarrays  except  OlA  have  been  affected  by 
occasional  communication  outages.  Problems  encountered  have  connection 
with  bad  communication  cables  in  the  ground  and  on  poles,  lack  of  power 
and  carrier  system  outages.  Also  other  irregularities  have  contributed 
to  errors  and  shorter  outages,  such  as  failing  CTV  equipment  like 
modems  and  SLEMS . 

In  addition,  there  were  problems  with  02B  (telemetry)  antennas  caused 
by  heavy  snowfall.  The  antennas  were  repaired  in  April. 

April  (weeks  9-13),  4.4  -  1.5.88 

02B  was  still  out  of  operation,  although  NTA/Hnmar  declared  the  line 
operational  13  April.  Danger  of  snowslides  prevented  the  power  plant 
people  from  repairing  the  broken  power  line. 

02B  (telemetry)  is  not  dependent  upon  external  power  (supplied  by 
solar  cells),  and  resumed  operation  (partly)  with  channels  23,  24,  2i 
and  27  active.  Channel  26  and  28  were  "dead”  but  the  NMC  staff  repaired 
the  antennas  (14,  15  and  25  April)  related  to  these  channels.  The 
antennas  had  been  broken  after  heavy  snowfalls  in  the  area. 

06C  was  affected  by  a  power  outage  week  15, 

The  performance  of  the  remaining  NORSAR  communications  systems  was  inost 
satisfactory. 


May  (weeks  18-22),  2.5  -  5.6.88 

02B  remained  down  also  in  May,  and  06C  had  its  modem/SLEM  checked  anti 
tested . 


The  remaining  communications  system  was  almost  free  from  trrors. 


June  (weeks  23-26),  6.6  -  3.7.88 

02B  resumed  operation  week  25  after  having  been  down  since  November 
1987. 


All  communications  systems  were  affected  for  about  1  hour  20  June  in 
connection  with  a  broken  coaxial  cable  at  Kjeller,  which  among  other 
things  carried  the  NORSAR  subarrays  and  NORESS.  OlA  and  OlB  resumed 
operation  after  a  couple  of  hours,  while  02C-04C  were  operational  a  tew 
houi's  later  in  the  afternoon,  06C  the  next  morning. 

On  15  June  the  IBM  2701  communications  adapter  failed.  This  machine 
interfaces  the  Modcomp  to  the  IBM  4341  via  two  highspeed  modems.  On  20 
June  the  source  of  the  failure  was  identified  and  the  problem  cor¬ 
rected.  In  the  meantime,  the  Modcomp  had  to  be  stopped/started  each 
rime  a  card  in  the  2701  was  replaced,  or  other  actions  implemented  as 
part  of  the  fault-finding  procedure. 

Apart  from  the  interruptions  mentioned  above,  the  commvinicat ions 
tv.stoiiis  were  reliable  in  June. 


July  (weeks  27-30),  4  -  31.7.88 

•Also  throughout  ttiis  period  the  systems  jul'T formance  vjas  most  satisfac 
tory.  An  exception  was  02B  where  NTA/Hamar  found  the  line  toward  the 
CTV  had  failed  and  06C  which  had  line  problems  in  connection  with 
1 ightning . 


August  (weeks  31-35),  1.8  -  4.9.88 

The  NTA  coaxial  carrier  cable  was  broken  19  August  causing  signiiicant 
problems  weeks  33  and  34.  All  subarrays  were  aflected  (OlA  for  a  very 
short  period) . 

01  p.,  which  had  been  down  since  15  July  due  to  bad  cable  was  back  in 
operation  1  September. 


The  remaining  systems  (except  for  02B,  03C)  were  back  in  operation  22 
August . 

02B,  03c  resumed  operation  later  due  to  a  SLEM-failure  (02B)  and  a 
cable  problem  (03C) . 


September  (weeks  36-39),  5.9  -  2.10.88 

Generally  satisfactory  performance,  in  spite  of  problems  such  as  a 
damaged  cable  after  lightning  (02B)  between  30  September  and  4  October, 
carrier  system  outages  affecting  03C,  04C  25  September,  and  between  30 

September  and  1  October. 

06C  was  affected  totally  85.71  minutes  weeks  37  and  38. 


O.A.  Hansen 


.APR 

a  V  ,  A  .  it  - 1  , 

MAY  (b) 

5)  (2. 5-5. 6) 

JUN  (4) 

(6. 6-3. 7) 

JUL  (4) 
(4-31.7) 

AUG  (5) 

(1  .8-5.9) 

SEP  (4) 
(5.9-2.10) 

AVERAGE 
1/2  YEAR 

0.02 

0.003 

0. 16 

0.02 

0.01 

0 . 04 

0 . 04 

d.02 

0.001 

0.15 

*58.93 

*89.64 

0.03 

0.05 

■■■'iuo.oo 

4'.-100 . 00 

*59,83 

0.03 

2)  0.02 

0.05 

0.03 

1' .  02 

0,002 

0.15 

0.02 

0.45 

0.03 

0.11 

0 . 02 

0 .  On  5 

0 .  1  7 

0.02 

0.01 

*1  .87 

0.35 

0.02 

0 . 09 

0. 16 

0.02 

0.01 

*1 . 83 

0,  35 

■  ^  '  0.001 

n  .  S  5 

2.6! 

*10.28 

*2  5.  2 

0 . 1  3 

0,97 

:P  14.29 

14.10 

9.03 

9 . 90 

16.85 

0 . 6 1 

0 , 2  / 

:.s  v)2B 

02F. 

02  B 

01R,O6C 

01R,06C 

02B,03C 

04  0 

01B,02B.n6i; 

0.02 

0.  14 

0 . 56 

0.02 

0.10 

0  .  1  3 

0.21 

St‘c-  s.  ft  ion  11.2  rogardiiig  h  i  giiros  profcdod  bv  .m  astoi'isk. 

V'ig'ii'c.s  represent; i iig  error  rate  (in  per  cent)  pre(.:eclt(i  liv  .i  nuniber 


c  .  ,  art 

'  related  to 

legend  i)»  low. 

rage  3 

week.s  , 

week 

14  N'/A 

Aver;ir.(> 

1  week.s, 

" .  > 

weeks . 

week 

33.34  N/A 

4  montlis 

4 

weeks , 

we'  k 

3  1  N/A 

8  ) 

3  niontli:, 

"  -  3 

weeks , 

week 

33,34  N/A 

9)  ... 

4  nionths 

"-  3 

weeks  , 

week 

?'<  N/A 

I'il*  Le  1 1  .  ?-  ■  1  Commun  i  ca  t  i  on.s  pe  r  toriu.ince  .  The  nnnibers  represent  error  ratr.'-;  i  i 
e.  .  t  ba.'.ed  on  ’otal  t  ran.siiii  t  ted  fr.aiiie.e /veek  (1  April  -  10  September 


11 


11,3 _ Event  Deteccion  operation 

In  Table  II. 3.1  some  monthly  statistics  of  the  Detection  and  Event 
Processor  operation  are  given.  The  table  lists  the  total  number  of 
detections  (DPX)  triggered  by  the  on-line  detector,  the  total  number  of 
detections  processed  by  the  automatic  event  processor  (EPX)  and  the 
total  number  of  events  accepted  after  analyst  review  (Teleseismic 
phases,  core  phases  and  total). 


Total 

DPX 

Total 

EPX 

Accepted  events 
P-phases  Core  Phases 

Sum 

Daily 

APR  88 

10825 

1135 

277 

72 

399 

11.6 

MAY  88 

5400 

941 

264 

96 

360 

11.6 

lUN  88 

f)600 

910 

199 

55 

2  ?  '♦ 

8.5 

JUL  88 

8850 

1191 

295 

8  7 

382 

12  ,  3 

AUG  88 

12750 

1191 

2/0 

71 

3^4  1 

11.0 

SEP  88 

9075 

945 

222 

60 

282 

9  4 

1527 

44l 

1963 

10,7 

Table  II . 3 . 1  Detection  and  Event  Processor  statistics,  October  1  987  - 
March  1988. 


B.  Paul. sen 


Ill . 


OPERATION  OF  NORRSS  AND  ARCESS 


1 1 1  . 1 _ Satellite  t-ransiiiission  of  NORESS  data  to  the  U.S. 

The  satellite  transmission  of  data  to  the  U.S.  from  the  NORESS  field 
installation  lias  fencraly  been  very  stable,  exept  for  occasional 
interruptions  due  to  power  break.s  and  control  line  breaks.  Tlitse 
out.ige  periods  ire  listed  in  Table  III. 1.1. 


10 

.Apr 

0600 

to 

0745 

power 

break 

20 

Apr 

0605 

to 

1310 

powei* 

break 

22 

Apr 

0617 

to 

0916 

powe  r 

l>reak 

;0 

Apr 

08  38 

1417 

poix'fc-r 

break 

May 

2056 

's  O 

5  .'•lay 

0023 

powe  L* 

break 

May 

2100 

t  ( ' 

10  May 

03.11 

power 

break 

1 1 

Hay 

220CI 

Vo 

2317 

power 

break 

.'O 

May 

1123 

to 

1124 

control  lines 

down 

2  0 

.1  un 

1050 

to 

112  3 

cent  rol  1 i nes 

down 

*  .1 

Aug 

1200 

to 

1205 

check 

ing  out 

transmi t  ter 

III. 

1,1, 

tage  pt 
1  ''88. 

■  r  i  0(1 

f('i  \(if 

'.ate  1  1  i  t  e  :  ransm 

■  pril 

-  S. 

'pten.be  r 

'he  total 

uptinu  for  the 

N''RESS  i:.Tr;h 

St  at  ioti 

tor  sa  t  e 1 11 

Sion  of  rlata  to  tlie  U.S.  '..’as  90.1*  .i;;  <ompared  to  97,/*  for  the 
previous  p-'-riod. 


.2 


Recording  of  NORESS  data  at  NDPC.  Kiellcr 

As  can  be  seen  from  Table  III. 2.1  the  reasons  causing  NORESS  outage  can 
be  placed  under  the  following  four  groups;  Transmission  line  failure, 
power  failure  at  HUB,  power  failure  at  NDPC  and  hardware  maintenance  or 
failure . 

The'  average  recording  time  was  97.8%  as  compared  to  95.4%  for  the 
previous  period. 


Date  Time  Duration  Cause 


18 

Apr 

1403-1427 

24 

m 

20 

Apr 

0817-1115 

4 

h 

58 

m 

22 

Apr 

0826-1047 

2 

h 

21 

m 

JO 

Apr 

IOdO-1736 

6 

h 

46 

m 

4 

May 

2311- 

49 

m 

5 

May 

-0458 

4 

h 

58 

m 

9 

May 

2313- 

47 

m 

10 

May 

-0509 

5 

h 

9 

m 

10 

May 

0722-0749 

27 

m 

11 

May 

0844-1304 

4 

h 

20 

m 

11 

May 

1308-1309 

1 

in 

11 

May 

1417-1535 

1 

h 

18 

m 

29 

May 

1626-1846 

2 

h 

20 

m 

29 

May 

1856-1911 

15 

m 

10 

Jun 

0501-0615 

1 

h 

14 

m 

20 

Jun 

1008-1908 

9 

h 

0 

m 

22 

Jun 

1336-1346 

10 

m 

28 

Jun 

1900-1932 

32 

m 

29 

Jun 

0825-0957 

1 

h 

32 

m 

1 

Jul 

0156-0304 

1 

h 

8 

in 

6 

Jul 

1016-1254 

2 

h 

38 

in 

14 

Jul 

2002-2106 

1 

h 

4 

m 

15 

Jul 

1810-2017 

2 

h 

/ 

in 

16 

Jul 

1658-1813 

1 

h 

15 

m 

19 

Jul 

1653-2037 

3 

h 

44 

m 

5 

Aug 

0248-0404 

1 

h 

16 

m 

10 

Aug 

1046-1102 

16 

m 

16 

Aug 

1330- 

10 

h 

30 

in 

17 

Aug 

-1102 

11 

h 

2 

in 

Transmission  line  failure 
Power  failure  at  HUB 
Power  failure  at  HUB 
Power  failure  at  HUB 

Transmission  line  failuie 
Transmission  line  failure 
Power  failure  at  HUB 
Power  failure  at  HUB 
Power  failure  at  HUB 
Hardware  maintenance  at  NDPC 
Hardware  maintenance  at  ^;DPC 
Hardware  maintenance  at  NDPC 
Power  failure  at  NDPC 
Power  failure  at  NDPC 

Transmission  line  failure 
Transmission  line  failure 
Transmission  line  failure 
Transmission  line  failure 
Hardware  error  at  NDPC 

Power  failure  at  NDPC 
Hardware  maintiuiance  at  NDPC 
Transmission  line  failure- 
Power  failure  at  NDPC 
Transmission  line  failure 
Power  failure  at  NDPC 

Power  failure  at  NDPC 
Hardware  maintenance  at  NDPC 
Hardware  failure  at  NDPC 
Hardware  failure  at  NDPC 


13 

Sep 

1630-1842 

2. 

h 

12 

m 

Transmission  line  failure 

13 

Sep 

1845-1846 

1 

m 

Transmission  line  failure 

15 

Sep 

1054-1251 

1 

h 

57 

m 

Power  failure  at  NDPC 

20 

Sep 

0903-1239 

3 

h 

36 

m 

Hardware  maintenance 

at 

NDPC 

2G 

Sep 

2303- 

57 

m 

Transmission  line  failure 

27 

Sep 

-0554 

5 

fi 

54 

in 

Transmission  line  failure 

27 

Sep 

1402-1406 

4 

nt 

Hardware  maintenance 

at 

HUB 

27 

Sep 

1445-1451 

6 

m 

Hardware  maintenance 

at 

HUB 

29 

Sep 

1124-1146 

22 

m 

Hardware  maintenance 

at 

HUB 

29 

Sep 

1203-1238 

35 

m 

Hardware  maintenance 

at 

HUB 

29 

Sep 

1324-1329 

5 

in 

Hardware  maintenance 

at 

HUB 

2  9 

Sep 

1409-1410 

1 

m 

Hardware  maintenance 

at 

HUB 

29 

Sep 

1415-1416 

1 

m 

Hardware  maintenance 

at 

HUB 

29 

Sep 

1421  -  1422 

1 

m 

Hardware  maintenance 

a  t 

HUB 

'"ah  1  0  I T  T  ■  2  .  I .  Int.f  rr  upt  i  on;;  in  NORRSS  recordings  at  NDPC,  April 
-  .September  1988. 


.'ionti'.iy  uptiiiio.s  f.'or  tiie  Noress  on-line  data  recording  task,  t, iking 
into  acc.onnt  all  tactors  (field  installations,  t ransniiss ion.s  line, 
da*,,  cent.,  r  operation)  affecting  this  task  were  as  follows: 


April 

98.0% 

day  : 

97 . 3% 

June  : 

00 

’’.ly  : 

98  .  * 

August 

96 . 9% 

September 

97.8% 

Fig.  III. 2.1  shows  tlie  upt  i/ne  for  the  i.lat.i  lerording  ta.sk,  oi' 
■'alently,  the  av.a  i  1  al)  i  I  i  ty  of  NORE.SS  data  in  our  tape  archive, 
d.'ty  - 1  iv  -  d.a  V  Vjasis,  for  the  reporlinp,  period. 


on  a 


15 


O  <  2  5  *»  5  Ji  T  8  H104  ' 


rn?«T2Q2r  i2.'5l’«2^2S?U‘'  726C'^r»  i 

UAW^ 


lAgu — IJI-i  2.1.  NORESS  data  record  i  nj'^  upt  imt'  tor  April  (top),  H; 
(middle)  and  June  1988  (boitnin). 


e>EP:  CEN'  CEN" 


^  T- 


£  S 


^  S  t  IB  9lO»  ^1215^^1S16t■n81a!>Q?1^2^^'^5^R2^^S^9R05 


Fig.  Ill  .2.1.  (cont.)  NORESS  recording  uptime  Cor  July  (top) 

August  (middle)  and  September  I  "'RR  (bottom). 


III. 3 


Recording  of  ARCESS  data  at  NDPC.  Kieller 


Monthly  uptimes  for  the  ARCESS  on-line  data  recording  at  NDPC,  taking 
into  account  all  factors  (field  installations,  transmissions  line, 
data  center  operation)  affecting  this  task  were  as  follows : 


April  : 

97.2% 

May 

79.6% 

June 

19.6% 

July  : 

87.8% 

August 

87.0% 

September 

92.0% 

The  main  reason  causing  most  of  the  ARCESS  outage  in  June  was  a  serious 
breakdown  of  the  Sun  2  processing  system  and  subsequent  work  in 
aonnection  with  transferring  the  data  recording  and  processing  to  a  Sun 
3-based  system.  In  August  work  on  a  new  power  line  at  tht  array  site 
was  the  main  reason  for  the  outages.  Other  reasons  causing  outage  in 
the  period  are:  Transmission  line  failure  or  line  testing,  power 
failure  at  the  HUB  or  at  NDPC,  hardware  and  software  work  or  tests  at 
NDPC. 

The  average  recording  time  for  the  period  was  7/. 2% 

Fig.  III. 3.1  shows  the  uptime  for  the  data  recording  task,  or  equi¬ 
valently,  the  availability  of  ARCESS  data  in  our  tape  archive,  on  a 
day-by-day  basis,  for  the  reporting  period. 


1  .  T<-i  st  ve  i ; 


IV. 


IMPROVEMENTS  AND  MODIFICATIONS 


IV  .  1 _ NORSAR  Detection  Drocc.s.«;inR 

The  'new'  NORSAR  detection  piocessor  has  been  running  satisfactorily 
during  this  reporting  period.  Detection  and  data  quality  reports  are 
now  grouped  into  day-files.  File  names  are,  e.g.,  NA088231.DPX  for 
detections  and  NA088231.REP  for  data  quality  reports  for  day-of-year 
231  1988.  File  naming  is  automatic  with  file  name  prefix  equal  to 
array  code  (i.e.,  NAO  for  NORSAR,  NRS  for  NORESS ,  FRS  for  ARCESS,...) 

Detection  processing  for  NORSAR,  NORESS,  ARCESS,  FINESA,  GRAFENBERG 
and  other  data  sources  has  been  executed  satisfactorily  with  the  new 
detector  program  using  a  macro  language  recipe  as  described  in 
scfjii  annual  report  No.  2-86/8  7. 


i  _ .MODCOMP  subarrav  communication 

No  irodi f ication  has  been  made  to  the  MODCOMP  system. 


I V .  3 _ NORSAR  event  process inp. 

Th.f  re  are  no  changes  in  the  NORSAR  evt-nt  lu'oct  .ssor  code. 


I V  ■  9 _ NORESS  detection  processing 

There  are  no  changes  in  tiie  'online'  RONAPP  detection  processing. 
Parallel  detection  processing  with  various  beam  sets  has  been  perforit 
to  itivesligate  detection  capabilities. 


IV  .5 


ARCESS  detection  processing 


Detection  processing  and  f-k  analysis  for  each  detection  have  been 
performed  in  an  off-line  mode  for  the  periods;  1987  (days  306  through 
364);  1988  (days  006  through  131  and  days  229  through  233).  Regular 
detection  processing  in  near  real  time  has  been  performed  on  the  Sun 
system  since  day  223  of  1988.  The  automatic  data  quality  control 
software  has  been  updated  to  detect  and  mask  channels  with  data  spil:e.s 
and  events  within  or  very  close  to  the  ARCESS  array. 


IV  .  6 _ Event  processinp 

A  new  event  processor  code  is  under  development  .  and  f-k  analysis 
for  NORESS/ARCESS/FINESA  has  been  periodically  performed  on  I'oth 
IBM  and  Sun  systems. 

The  event  processor  package  follows  the  same  design  as  the  detection 
processor  with  regard  to  macro  language  inptit .  The  ARRMAN  program 
system  which  handles  input  from  NORSAR,  NORESS ,  ARCESS,  FINESA, 
Grafenberg,  GSE  level  2  as  well  as  the  CSS  2.8  format,  is  included  'n 
both  packages.  The  package  operates  on  archive  tapis,  disk  files  and 
disk  loops,  as  available.  CSS  2.8  is  implemented  only  on  the  Sun 
computers,  whereas  the  other  formats  can  be  accessed  from  eithi'i"  IBM  or 
Sun . 

NOGRA  -  NORSAR  Graphics  package  -  has  been  implemented  by  NTNF/bORSAR 
for  graphics  output  under  IBM/GDDM,  IBM/GKS,  IBM/PIilGS,  Sun/Xll  and 
Sun/Laser.  Standard  graphics  user  routines  using  NOGR.'i  have  been 
developed  for  f-k  plots,  power  spectra,  interactive  trace  display.s, 
etc . 

The  event  processor  code  executes  broadband  f-k  analysis  and  uses 
velocity  to  report  preliminary  phase  determination.  The  intention  is  :o 
process  detection  information  from  one  array  at  a  time,  and  produce 
reports  with  velocity,  azimuth,  quality,  polarization,  frequency. 
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amplitude,  and  (preliminary)  phase  determination.  Preliminary  epicenter 
determination  will  be  based  upon  one-array  data.  Improved  location 
estimates  will  be  based  upon  these  reports  using  several  arrays  and 
Scandinavian  stations  when  available. 

Emphasis  will  be  put  on  flexibility  and  suitability  for  research 
applications,  and  the  work  will  be  coordinated  with  the  efforts  on  the 
Intelligent  Array  Processing  System  (IAS),  which  is  scheduled  for 
delivery  to  NORSAR  during  the  summer  of  1989. 


IV .  7 _ Upgrade  of  the  ARCESS  data  acquisition  and  processing 

hardware 


I'p  to  June  1988,  the.  ARCESS  data  acquisition  system  at  NDPC  was  based 
on  a  temporary  solution,  using  a  Sun-2  computer  on  loan  from  Science 
llorir.ons.  Vital  hardware  components  of  this  system  failed  in  early 
June,  and  there  was  a  halt  in  the  recording  of  ARCESS  data,  which 
la.sted  until  the  new  Sun-b,ased  system  (wViich  was  delivered  in  late 
June)  became  operational  in  early  July.  After  the  installation  of  the 
ncv.'  system,  the  hardware  configuration  Is  essentially  as  shown  in 
Fig.  IV. 7.1  in  NORSAR  Scientific  Report  No. 1-87/88.  Two  Sun-3/280- 
corrputers  handle  the  data  acquisition  and  data  processing  tasks,  and  a 
Sun-  3/260-compu:er  is  available  for  off-line  data  analy.sis. 

J .  Fyen 

R.  Paulsen 


V. 


MAINTENANCE  ACTIVITIES 


V. 1 _ Activities  In  the  field  and  at  the  maintenance  center 

This  section  sununarlzes  the  maintenance  activities  In  the  field,  at  the 
maintenance  center  (NMC)  at  Hamar  and  NDPC  activities  related  to 
monitoring  and  control  of  the  NORSAR,  NORESS  and  ARCESS  arrays. 

Activities  related  to  the  NORSAR  array  have  been  diverse,  and  most 
tasks  concern  corrective  maintenance.  Scanning  Table  V.l,  we  find  that 
repair/splicing  of  seismic  cables  is  prominent,  but  we  also  find 
repair/replacement  of  electronic  equipment  such  as  modems,  SLEMs  and  LP 
seismometers  adjustment  devices  (RCDs) . 

In  addition  to  repair/repositioning  of  antennas  related  to  the  02B 
telemetry  stations,  the  NMC  staff  have  taken  care  of  such  jobs  as 
adjustment  of  Long  Period  seismometer  parameters  Free  Period  (FP)  and 
Mass  Position  (MP) . 

The  NDPC  staff  regularly  monitor  the  subarray  electronics,  including 
communication  systems.  LP  seismometer  parameters  FP  and  MP  are  remotely 
adjusted.  SP  channel  parameters  such  as  RA-5  gain,  RA-5  3  dB  point, 
filter  ripple,  LTA  time  constant,  seismometer  .sensitivity  and  natural 
frequency  are  evaluated  by  means  of  special  on-/of f  1- i''G  programs. 

The  NORESS  field  system  performed  entirely  satisfactorily  throughout 
the  reporting  period,  and  only  a  few  corrective  actions  were  needed. 
Details  are  given  in  Table  V.l. 

During  a  visit  9-15  June  to  ARCESS,  a  new  Global  Positioning  System 
(GPS)  Synchronized  Clock  was  installed  by  Sandia  engineers  in  coopera¬ 
tion  with  NORSAR  field  personnel.  During  the  same  visit,  all  the  fiber 
optical  cards  were  modified.  The  modification  made  it  possible  to 
adjust  the  frequency  and  the  phase  of  the  clock  signal  for  the  optical 
data  links  between  the  HUB  and  the  remote  sites. 


The  KS- 36000  borehole  seismometer  has  not  been  operating  properly  due 
to  120  Hz  oscillations  in  the  two  horizontal  sensor  modules.  It  was 
therefore  decided  to  pull  the  seismometer  from  the  borehole  and  to  send 
it  back  to  Teledyne  Geotech  for  repair.  This  was  done  during  the  9-15 
June  visit  to  ARCESS. 

During  the  same  visit,  a  cold  soldering  on  the  interface  card  in  the 
Communications  Interface  Module  (CIM)  was  repaired.  After  that  repair 
was  carried  out,  there  has  been  no  failure  in  the  data  link  to  Kjeller, 
which  can  be  traced  to  the  CIM  at  ARCESS. 

It  was  found  in  September  that  the  GPS  unit  did  not  operate  properly. 

We  found  that  the  United  States  Naval  Observatory  (USNO)  and  Department 
of  Defence  (DOD)  had  changed  the  data  content  of  the  GPS  signals  being 
transmitted  by  satellites  6  and  9.  This  change  of  data  contents 
inhibits  the  NAVCORE  Signal  Processor  from  acquiring  these  satellites. 
In  addition,  this  change  also  inhibits  the  NAVCORE  from  acquiring  any 
other  satellite  while  satellite  6  and  9  are  selected. 

Since  the  USNO  was  in  the  process  of  changing  the  data  content  on  all 
of  the  satellites,  and  in  order  to  fully  correct  this  problem,  the  GPS 
unit  was  sent  back  to  Kinemetrics  for  a  software  modification  that  will 
allow  the  NAVCORE  to  decode  the  new  data  format. 

In  order  to  have,  a  timing  system  at  ARCESS  during  the  period  the  CPS 
will  be  in  the  U.S.  for  repair,  the  two  I.F-DC  Synchronized  Digital 
Clocks  were  modified  to  work  properly  under  marginal  receiving 
condi tions , 

During  the  summer  months,  all  seismometer  housings  of  the  ARCESS  arrav 
were  covered  by  moss,  in  order  Co  reduce  as  much  as  possible  on 
background  seismic  noise  from  wind  and  rain. 


Subarray 

Area 


Task 


Date 


02B 

(telem. ) 

Temporary  antenna  repair  and  repositioning 

14,15,25 

April 

06C 

Three  visits  to  the  SA  in  connection  with 
testing  of  a  SLEM  which  had  caused  timing 
problems  before 

May 

NDPC 

Remotely  measured/adjusted  l.P  seismometers  with 
respect  to  Mass  Position  (MP)  and  Free  Period 
(FP) .  Subarrays  OlA,  OIR,  02C,  03C  and  04C 

May 

02C,  03c  SP  channels  evaluated  by  means  of  the 
Offline  program  Chanev  SP.  Parameters  such  as 
filter  ripple,  LTA  time  constant,  RA-5  gain, 
RA-5  lower  3  dB  point.  Seismometer  sensitivity, 
and  natural  frequency  determined 

02B 

Adjusted  SP/LP  gain  and  Offset 

3  June 

04C 

Remote  Centering  Device  (ROD)  EW  so-israome  t  er 
replaced 

24  June 

OLA.OIB 

SP/LP  gain  and  offset  adjusted 

27  June 

OIB 

Seismic  cable  SP02  spliced 

29  June 

03C 

SP/LP  gain  and  offset  adjusted 

3  0  ,1  une 

02P. 

Found  seismic  cable  SP02  damagtHt 

30  Ji.ire 

NOR ESS 

Array  checked  together  with  Sandia  repre.^tu- 
tatives 

8  Juno 

ARCESS 

NMC  staff  and  Sandia  represent.at  i  ves  visi'ed 
the  array  in  connection  with  modification  and 
timing  of  Fiber  Optical  systems.  A  Global 
Positioning  System  (GPS)  was  installed 

10-1.5  June 

Table  V.l, 
including 

Activities  in  the  field  and  the  NORSAR  inaintr 
NDPC  actitivites  related  to  the  NORSAR  array,  1 

nance  vente'. 
Apr i 1  - 

SO  September  1988. 


Subarray 
.■-Ire  a 

Task 

Date 

h'DPC 

Daily  check  of  SP  and  hP  data. 

including  comm.  June 

systems.  LP  seismometers  outside  .specifications 
adjusted,  Free  Period  (FP)  and  Mass  Pos .  (MP) . 
SP  channels  OIB,  02B  and  04G  analyzed  with 
respect  to  vital  parameters  by  means  of  the 
Offline  program  Chanev  SP 


OlA 

Cable  repair  SP03 

3,8,12 

July 

0 1  B 

Cabl,-  repair  SPU? 

1  July 

DIB 

Cable  repai'.'  SPOj 

19  July 

'■’211 

tiable  t'- pa  i  r  Sl'U/’ 

22  J  u  1  y 

(V*C 

Cable  r?  p.i  1 1:  Si'Oi 

20.26,27 
28  July 

OjC 

Cable  refiair  SP04 

13,14  Jul 

:;dpc' 

Daily  routine  che-'ks  throughout  the  month 

July 

ARCF.SS 

Replaced  satellite  receiver  clock  hv  the  "old" 
i  ecu*  i '.’er.-;  .t!F  1  array) 

4  .August 

N’ORl-'SS 

Satellite  t.  •'ansini  t.ter  frequency  decresased  by 

305  lit.  Hub  previ-ntive  nia Lnt enance  carried  out. 
Romoti  .site-,  power  supply  rep.iiiu'd 

19  .August 

24  .-August 

'!  i  B 

Cable  work  SPO.') 

3 . 4  Augu.s 

02  b 

Repj.'a'ftl  mo.lem  ana'  SI.FM  Digital  Unit 

11,17,2). 
24  .’vugust 

02P, 

Caljle  work  SPO? 

10  .Xni’.ust 

04  C 

The  S.\  Was  visited  in  ordtr  to  solve  a  .special 

NS  LP  seism,  proljlein 

19  August 

I'able  V.l  ( c  o  n  r  .  ) 


Subarray 

Ar€!a 


Task 


Date 


NDPC  Regular  check  of  SP/LP  data  and  comm,  systems. 

Weekly  calibration  oi  SP/LP  seism,  carried  out. 
Adjustment  of  LP  seismometer  (MP/FP)  done  when 
outside  tolerances.  In  August  the  capability  of 
the  02B,  02C,  03C  and  04C  A/D  converters  was 
verified  by  "online"  data  acquisition  "test  48" 
and  the  "offline"  program  MISNO 

NORESS  Routine  visits 

ARCESS  On  site  A1  the  Optical  Fiber  Transmitter  was 
replaced.  On  D6  the  Hub  61  and  70  cards  were 
replaced.  The  DHL  70  card  on  the  north -south 
component  on  C2  was  replaced. 

NORSAR  Regular  check  of  SP/LP  data  and  comm,  systems. 

Weekly  calibration  of  SP/LP  seismometers.  Mass 
Position  (MP)  and  Free  Period  (FP)  outside 
tolerances  adjusted.  Besides,  data  acquisition 
"test  48"  and  offline  data  analyzing  program 
Misno  run  on  subarrav  01 A 


Tab 1 e  V . 1 .  (cont.) 


V .  2 _ Array  status 

No  changes  or  modifications  have  been  implemented  since  ti' 
report . 

As  of  30  September  1488  the  following  NORSfXR  channel. s  de-  l 
tolerances : 


OlA 

01 

8  Hz  filter 

02 

8  Hz  filter 

04 

30  dB  attenuation 

UIB 

05 

biiti  cable 

August 


Set  tembc' r 
September 

September 


c  I;,.-:. 

,iLed  from 


0.?B  08  FP  noL  .idjiist  ablt'/nif ‘T 1 1: xn  Nl)!'  ' 

L'F  .'8  FP  uoP  .idj  iiy ill  1  >■  from  NDPi' 

Ob'"  1)8  liroadli.ind  I  liter  i  ll.‘■.l  al  i  ed 

On  the  oceaslon  of  flie  Semipalat  i  tisk  JVE  explosion  on  14  September,  Liu 
gain  of  the  NORESS  1-componoitt  instrument  at  site  C2  was  reduced  by 
dB .  This  change  was  effective  on  13  September  at  1700  GMT,  and  the  gain 
was  back  !o  normal  from  16  Sept;ember  at  1400  GMT. 

No  7RCESS  ch,inn(.l.s  d; ia  t..  (.1  from  their  standards  during  the  reporting 
pe  r ; od . 


O.A.  dansen 
r.W.  Larsen 
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VII.  SIIMM.AilV  0)  TECHNICAL  REPO.'IPS  /  PAPERS  PUBLISHED 


VII  .  I  Si^eccral  analysis  of  Shap.an  Rtvor  explosioti.s  recorded  nt  NORSAR 
and  NORESS 

A.s  .shown  hv  Run-dal  ai  d  tlokland  (1987),  NORSAR  recordings  of  Lg  and  P 
cod, a  can  providi  ■.'(  rv  .stable  c.s t  imal.r.';  of  the  iiungnitudes  of  underground 
nuclear  e.xpLosions  a’  the  Shagan  River  te.st  site.  These  data  thus  hold 
cons  ideralile  promise  in  relation  to  obtaining  accurate  yifild  estimate.-:; 
for  the  purpose  of  verifvir.;',  ■■  tiiroshold  ttist  ban  tre.ity. 

To  ir.vesrigate  the  obsorv.it  i  ona  1  basis  for  estimation  of  magnitudes 
ba.sed.  on  P-coda  and  Lg  me.isuretitents ,  we  have  calculated  power  spectra 
fre-m  NORS.sR  recn  rd  i  inv"  of  about  twenty  Shag.an  River  oxplo.sions. 

The  mean  specf,  ui'  of  the  NORSAR  short-period  cliannels  and  the  cnrve.s 
to  ,1  ;■  -sent  ing  p  !  u.s  ,  ir  i  iv.s  two  standard  deviations  are  estimated  for  noi-.e 
pi'c-eding  the  1-pii-ise.  for  A-roda  and  for  Lg.  The  time  windows  are 
eqtiil  to  those  .ipplied  by  Ringdal  tind  ItokL.and  (1  987)  for  estimating  RMS 
ba.si-d.  mag:  itudes  We  will  r.ot  go  into  detail  on  their  procedure,  but 
j'oit't  out  that  the  RM!  of  r.oise.  P-coda  .and  l.r.  v./ere  c.ilculated  from 
t.r.ac,.,s  b.Ti'dpass  ■  t  i  '  '  e  i  fd  ijetwoin  0  .  C  and  ;.(i  I':'.  Fig.  VI I  .  1 , 1  ( a  -  c  )  gi'a 
the  spectra  froii  the  .loint  7e  r  i  f  i  c.a  t: ;  on  Fxperimont  (viVF)  explosion  ir. 
tilt  Sliaga.a  Riv’cr  rcgicii  of  Septe-mber  I A .  19.°'^.  TI'  ;  I  -i-p  mi.f-,;- 

spectrum  across  NORSAR  as  well  as  euivcs  corvts|iond  i  ng  to  plus/niiiius 
two  standard  deviations.  ii  tin  r  ''ig.-  0 . 0  - 1 . 0  H."  ,  we  find  that  the 
variations  across  t!u  NORSAR  ari.a-.-  art-  of  comparable  ;;ize  for  bocli 
noise,  P  coda  and  I.g ,  and  the  saive  v.iriai  ion  cli.iract  eri  sties  also  app’.v 
to  the  other  events  i  I'vtist  i  gated  in  Ibis  stufiv. 

A  procedure  to  compens.i te  for  the-  tiackground  noise  .,.evel  torms  part  ol 
the  RHS  magnitude  measurements  This  procedure  can  be  simulated  on  tht 
power  spectra  bv  sub  1 1  ac  t  i  ng,  the  pre-P  noise  from  the  P-coda  or  Lg 
spic!  r.a.  in  Fig.  V 1  I  .  1  .  7  we  illustrate  the  hackgroun'i  noise  comptnisa- 
t  ion  b’  showing  the  sptetra  uf  flu  noise,  of  Lg  and  of  Lg  minus  noise. 


We  no(<!  that  iti  tliL.s  the  tff.ovl  of  tlii-  iioiste  roiiii/ensa  l  ion  in  kih.i'i  1 

In  Llie  fie.qnenoy  band  with  the  maximum  power  (around  0.8  Hz).  In  fac  . 
tlie  noise  compensation  only  becomes  important  for  the  Iov;er  magnitude 
events.  In  the  following,  refercnct's  will  be  made  to  noise-corrected  L.g 
spectra  only. 

Note  that  the  Lg  spectrum  of  Fig.  VI 1.1. 2  exhibits  a  peak  betwetu  0 
and  0.8  Hz.  From  the  NORSAR  short -period  response  given  in  Fig. 

VII. 1.3,  we  can  see  that  in  the  frequency  range  0.6  to  3.0  Hz  the 
amplification  is  varying  by  a  factor  of  ten.  If  the  dominant  frequency 
of  the  Lg  phase  were  to  vary  significantly  from  event  to  event,  the 
RMS-based  magnitude  estimates  could  be  influenced  bv  the  'arying, 
amplification.  To  investigate  this  problem  further,  Fig  VIT.l.ia  shov.'s 
expanded  Lg  spectra  of  eleven  events  from  the  soutlr./e.c  re  i  n  part  of  tlie 
Shagan  River  test  region.  Tlie  RMS  Lg  ma.gnitudes  range  from  3.6?  to 
6.19.  Although  there  is  a  trtmd  of  lower  domin.int  f  rc  |ui  ici.  s  jor  tiu 
highest  peaks,  the  actual  varitition  is  small.  In  Fi-..  Vli  l.ab  ■..>.■  slew 
similar  spectra  for  seven  tvents  from  the  nor theast ■ in  part  oi  the 
Shagan  River  test  region.  The  RMS  l.g  magnitudes  for  these  events  varv 
from  5.87  to  6.11. 

Ringdal  and  Hokland  (1987)  found  in  their  stuuv  ot  l.g  .n.c.  i’  cod, 
magnitudes  from  the  Shagan  River  test  site  that  there  was  a  sig¬ 
nificant  regional  anomaly  within  that  site.  In  the  noithi  .aste  rn  p.irt  . 
the  P-coda  and  ISC  magnitudes  were  consistently  low  comp,!  ;e;'  to  Lg , 
whereas  in  the  southwestern  part  tliey  were  consistently  h  i  r,h .  i\i 
invc'Stigaee  whether  this  anomaly  is  reflected  in  tie  I.r  ,  ec  :  w.’ 

have  plotted  in  Fig.  VII.  1.6  the  peak  I  ref|uenc  ies  lor  tlu  events  iir.s 
tigated  as  a  function  of  RMS  Lg  magnitude,  with  diilcr-u:  svmho!  t vj .  . 
for  the  two  subregions.  Although  there  are  a  couple  ot  out  liers,  hotl; 
the  events  from  the  NE  part  (crosses)  and  the  SE  part  (.tilled  rec¬ 
tangles)  follow  the  same  trend  of  lower  dominant  f  iw  quenc ■ c s  for  higiur 
magnitudes.  From  the  results  given  in  Fig.  VII. 1.5  wo  can  infi ;  that 
the  Lg  peak  frequency  characteristics  do  not  differ  s i gni f Leant ly  fmii 
the  NE  to  the  SU  part  of  th»  Shagan  River  tost  site.  The  overall 
variation  of  the  dominant  frequency  among  th<'  ovt-nts  of  differeor 


inagpitiide'  is  less  rh..i!  0.':)  II2 .  l-s-  lliore  fore  liiid  Chnr.  .Tirp  1  i  f  Leaf  i  on 
di  Ef  eronca- s  ,  sc  ^  Ki;-',.  Vn.1.3,  '-.'ill  not  infliic’nof  t!io  va  r  i  .1 1  i  on  of  t!. 

iragiii  tudo  estiiPal.es  s i  gnif  leant;  1  v  . 

In  rig.  Vll.l.ta  and  Vll.l.rh  we  have  calc'ulat.ed  I'-codca  and  Lg  spectra 
for  two  Shagan  river  events  with  comparable  RMS  Lg  magnitudes.  The 
everiL  represented  in  Klg.  VI  I. 1.6a  is  located  in  the  SW  region  and  has 
an  roiS  Lg  magnitude  t'f  S.'hs,  whereas  the  event  given  in  Fig.  VII. 1.6b 
is  located  in  tlic  NE  region  and  has  a  magnitude  of  6.8/.  The  feature  we 
want  to  emphasi;;e  from  thi.se  figures  is  that  the  P-coda  spectrum  of 
the  SW  evc'nt  is  well  ihove  that  from  the  NE  c-vent  in  the  fntire 
frequenev  range  0.6  to  3.0  !iz .  On  the  other  liand,  the  differences  in 
the  Lg  spc'otra  .ire  sir.ill  and  ar.-  confined  to  the  frequency  1  ange  0.6  to 
L.O  Hr. 

For  events  lulo,.  .1  certain  magnitude,  the  SNR  of  i.g  at  NORSAR  i.s  too 
low  for  applicpition  ol  Lg, -based  magnitude  in  asurements.  In  theory,  the 
SNH  eiiuld  'C  in'nro”ed  through  lu-amformlng .  hut  in  practice  the  Lg  ph.ise 
ii.ns  t'lo  I',.-,-  coh<  reni,".'  .icros.s  NORSAR  for  tin's  to  be  meaningful.  On  th«‘ 
ot. }'.<■;■  he,;:d.  the  N'OKK.Ss  arriv  h/n;  .shown  an  excellent  captibility  of 
i:i  ;  .'ir.,'  'he  SNR,  In  Fig.  '.'I  I  .  1  .  /a  we  show  tlie  ine;ui  NORESS  spect  ra 

for  •'.oi^if  'preci  : : 'I  •,  fl,(.  Pphase  and  for  L;'. .  The  even:  considered  is  the 
JVE  .xpln;  inn  o;  Si: 'p :  •  nhe  1  ii,  .A:  O.S  lit.  tl.c  S'.'R  i ahou'  11 

decibels.  hv  foniuni’,  a  be.iir  from  the  center  instrument  ,'ind  the  D-rin,; 
of  tn.e  N'tid'.S.S  a;  ray  with  ,s  t ‘-e  r  i  tig  di-l.'ivs  eor  res|)ond  i  tig,  to  an  tip  pa  r.  ■  11 : 
vclocitv  of'  <1,3  l:m/',  .ind  the  tirimuth  to  the  Shag.an  River  test  site  180 
degree,-;),  the  Sl.'l:  e.aii  he  significant  ! '/  improved.  in  1' i  g .  VI  1  .  1  .  7ti  ,  w, 
slio',..  rile  h.-.ini  spectra  of  noise  ,ind  Lg,  .ind  find  and  SNR  oi  17  deeili'  Is 
/it  1.8  li;;  .  impl’.ing,  that  Ni'RESS  I^g  measurement.'-;  of  Shaggui  River 
e.xplo.sions  mai/  ire  deiie  for  ev«'nt  of  6  dB  (0.3  mg  i.inits)  lower  tluin  tie- 
NORSAR.  l!ov;ever,  in  the  low- SNR  e.a.se.s  where  we  have  to  apply  the 
Ireainf orniing  tecl  tti'iue  on  NORESS  array  recordings  to  obtain  an  Lg 
m/'gt:itndo  e  s  t.  i  ir/i  t  e  ,  ll'r  varianee  in  the  e.stimate  will  he  larger  than 
for  t Ik  high  SNR  cases  win  re  we  can  a"erag,o  over  the  full  NORSAR  arrav. 


To  illustrate  how  NORESS  bearaformliig  works  for  the  Lp,  phase,  v't-  have 
applied  the  wide-band  slowness  estimation  techinque  to  NORESS  short - 
period  recordings  of  a  16  min  long  wavetrain  comprising  all  phases  from 
the  JVE  explosion.  The  center  instrument  and  the  C-  and  D-rings  were 
analyzed  in  the  frequency  band  0.6  to  3.0  llz.  Each  time  window  was 
three  seconds  long  and  the  separation  between  the  windows  was  one 
second.  The  results  are  given  in  Fig.  VII.  1.8  and  show  the  folltiwing 
In  the  upper  panel,  the  intermediate  period  vertical  ch.e.nel,  bandpa.s^: 
filtered  between  0.6  and  3.0  Hz  is  displayed.  In  orcier  ;  inoro  clearly 
visualize  the  PP  phase,  occurring  after  iibout  120  seconds  and  the  Lg 
phase  arriving  between  750  and  870  seconds,  we  havt  clipped  the 
amplitude  of  the  P-phase  in  the  plot.  In  the  second  panel,  th  aziiir.ith 
from  the  slowness  analysis  is  given.  The  size  of  tlie  ri'.'ps  repiesen;  a 
coherency  measure  of  the  slowness  solution.  Although  there  is  a 
relatively  large  scatter  in  azimuth  around  the.  theoretical  value  of 
about  80  degrees,  the  time  interv.tls  arountl  P.  the  enrlv  P-cc>da,  PP  ar.d 
T.g  show  a  more  uniform  pattern  than  the  r^'sL.  The  slowness  or  anpnr<;.‘ 
velocity  estimates,  given  in  the  lower  panel,  show  that  hetwcei:  P  an’ 

PP  the  apparent  velocity  is  consisteittly  above  10  km/s .  It  ihn  droi  s 
to  ;ielow  5  km/s  after  about  6  minutes,  and  then  stays  tit  about  rhi.s 
level  throughout  the  wavetrain.  The  relatively  conslstei.r  azin;  ith  ti  id 
apparent  velocity  estimates  within  the  Lg  wavetrain  expltiin  whv  the 
beamforming  works  well  for  this  phase. 

Conclusions 

In  this  study  we  hav<^  prt'sonted  some  results  i  1  ]  usti'.at  i  nr  ..oiiu  ■)  thi 
features  related  to  magnitude  estimation  based  on  RIIS  l.r  uirl  RI-LS 
P-coda  measurements.  The  Lg  spectra  from  both  the  NL  and  'la  ,S\'  [i.iri  of 
the  Shagan  River  region  show  little  variation  in  si''Ctra;  vliape  and 
dominant  frequency.  Even  though  the  RMS  Lg  inag,nitudes  .,t  <  computed 
from  traces  filtered  in  the  0. 6-3.0  Hz  band,  the  spoc'ra  uliow  iliat  tla: 
signal  energy  level  between  0.6  and  1.0  Hz  essentially  det  oiinin.-s  the 
Lg  magnitudes. 


For  event;',  from  ''I'v'  Sluigau,  the  .spectral  difference  between  P-coda  and 
Lg  is  larger  than  for  events  from  the  NE  part.  Thi.s  applies  to  the 
entire  frequenc'/  rarige  0.6  to  3 . 0  Hz  .  It  also  follows  that  the 
spectral  level  of  the  entire  frequency  band  analyzed  (0.6  to  3.0  Hz) 
contr  ibtites  to  t  lie  RMS  P-c.oda  magnitude  estimates. 

In  cases  where  the  single  station  SNR  of  the,  bg  is  "oo  low  for  Lg 
magnitude  estimation,  we  can  employ  the  beamforming  capability  of  the 
NORESS  array  to  improve  the  SNR.  About  6  decibels  SNR  improvement  can 
be  achieved,  i.e.,  aliout  0.3  magnitude  units.  The  applicability  of 
beamforiring  of  the  I.g  ph.ises  has  been  demonstated  by  running  moving, 
time  window  slo'*’uc;ss  analysis  on  a  16  minutes  long  window  cov'ering  all 
phases  from  a  Shagan  River  nuclear  explosion. 

T.  Kv*rna 
F.  Ringdal 
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Ftp.  VII ■ ] ■ 1 .  Mean  uncorrecred  NORSAR  power  spectra  of  noise  precedir, 
the  P-phase  (a),  P-coda  (b)  and  Lg  (c).  The  upper  and  tower  curvets 
indicate  plus/minus  two  standard  deviations  The  event  .inalyzed  is  th 
Joint  Verification  Experiment  (JVE)  explosif)ii  at  the  Sltagan  River  test 
site  on  September  lA,  1988. 
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P~ i  g  .  V 1 1 .  1 . 2  ■  The  noise  compensation  procedure  is  illustrated  in  this 
figure.  The  upper  spectrum  represents  the  Lg  phase  from  the  JVE,  the 
ln;rcr  .spectrum  represents  noise  preceding  the  P-phase  and  the  dif¬ 
ference  is  given  in  the  middle.  Note  that  the  frequency  range  around 
the  spectral  peak  of  Lg  is  marginally  influenced  by  the  noise  compensa 
tioii . 
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'•'ifi  VJ  t  .  1 .  3  ■  NORSAR  short-period  velocity  response  function. 
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Fig.  VII. 1.4a  Expanded  plot  of  Lg  spectra  calculated  iroiii  NORS.^R 
recordings  of  eleven  explosions  from  the  southwestern  p.u  t  of  the 
Shagaa  river  test  site. 
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Fig.  VII. 1.4b.  Expanded  plot  of  Lg  spectra  calculated  from  NORSAR 
recordings  of  seven  explosions  from  the  northeastern  part  of  the 
Shagan  river  test  site. 
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Fig..  VII  ,1,5.  Peak  frequency  ver.sus  RMS  Lg  magnitude  for  the 
given  in  Fig.  VII. 1.4a  and  4b.  Crosses  represent  events  from 
part  of  the  Shagan  River  test  site,  wliereas  filled  rect.angles 
event,s  from  the  SW  part:. 
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Fig.  VII.l,6a.  P-coda  and  Lg  spectra  from  the  JVE  explosion  of 
September  14,  1988.  This  event  is  located  In  the  SW  part  of  the  Shagan 
River  test  site.  The  RMS  Lg  magnitude  is  estimated  at 
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Ilfc.  vn.l  6b.  p-coda  and  Lg  spectra  from  an  event  of  12  Dei  ember 
1984.  This  event  Is  located  in  the  NE  part  of  the  Shagan  River  test 
site  and  the  RMS  Lg  magnitude  is  estimated  at  5.87. 
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LL£j — YLl-r  1,7.  Illustration  of  SNR  gain  for  Lg  by  beamforming  using 

NORESS  data  from  the  JVE  explosion.  Panel  a)  shows  uncorrected  NORESS 
power  spectra  for  Lg  and  noise  averaged  over  all  individual  seis 
mometers.  Panel  b)  is  based  on  a  beam  formed  from  the  center  instrument 
and  the  D-ring,  using  steering  delays  typical  of  Lg  phases  from 
Semlpalatinsk  (phase  velocity  A. 3  kra/s,  azimuth  80  deg).  Note  the 
considerably  greater  SNR  for  the  NORESS  beam. 
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Fig,  VII. 1.8.  Results  from  slovmcss  analysis  of  NORESS  r<-cordings  of 
the  JVE  explosion.  The  upper  panel  shows  the  intermediate  period 
vertical  component  bandpass  filtered  between  0.6  and  3.0  Hz.  The 
middle  panel  gives  the  azimuth  solutions  and  the  lower  panel  the 
absolute  slowness  (inverse  of  apparent  velocity)  solutions.  The  size 
of  the  circles  represents  a  coherency  measure  of  the  indi- idual 
estimates . 


VII  .  > 


Statistics  of  ISC  r-r.ivel _ t iiric  rcsiciiials 


The  location  of  scisii.ic  events  hy  a  network  of  stations  requires 
adequate:  tiieoretical  travel  tiiiK-s  are  available,  either  bv  interpola- 
tioi;  in  tables,  or  by  calc.n  ating  the  times  in  a  reference’  velocity 
model.  Despite  obviou.s  shortcoming's  the  Jef  f  rey'' -  Bui  ion  tables  arr- 
stiil  in  use  at  the  major  seismological  centers  for  locating  events. 

The  more  recent  PRfM  piode  1  (Dziewoiiski  and  Anderson,  1981)  is  more 
satisfying  in  tiiat  i wa.s  constructed  to  fit  a  large  seismological  data 
set  including  free  oscillation  c igenfrcquencies ,  but  the  original 
transversely  isotropic  model  is  not  well  suited  for  routine  travel  time 
calculations,  and  the  isotropic  version  of  PREM  ha.s  not  been  adequately 
tested  against  .irriva'i  rime  observations.  Here  wc  report  on  the 
s  tat  i s 1 1  c  i  of  tt' 1  e.sc  i.smic  travel  time  residuals  with  respect  to  the 
L.sorropLc  PREM, 

Ve  li.ive  extract,  d  P  ai  rival  time  data  from  the  ISC  bulletins  for  the 

198A  ,  peP  and  PKP  for  the  years  1975-198a,  and  PKKP  for  the  years 
1964-198A.  Additional  PKKP  and  PnKP,  n  >  2,  were  taken  from  bulletins 
of  the  original  IWSA  and  NORSAR  arrays,  and  from  special  publications 
All  data  were  s  ;bjected  to  a  standard  processing  sequence,  similar  to 
that  of  others:  Residuals  were  computed  relative  to  PREM,  subjected  to 
station  corrections,  cori'ecttd  for  .^'l  1  IptLci  tv  and  lower  mantle 
variations,  and  corrcctccl  lor  the  cift'ct.':  of  soui'ce  structure  atui/or 
mis  1  oca t i on .  Data  belonging  to  a  particular  branch  were  finally 
averaged  to  form  'summary  ray'  data,  ba.sed  on  piiirs  of  approx  ima  te  ]  v 
equal  area  blocks  (equalling  10  >;  lo"  ;,t  t  iu-  equator).  For  detail.s  of 
the  data  selection  and  processing  we  rcfi  r  to  Doornb.r.s  and  Hilton 
(1988).  The  numljer  of  'summary  rav'  data  ILnallv  obtained  were,  .'^al'l  fe.r 
P,  lt68  for  PcP.  139')  for  PKP  (lie),  8,/]  for  PKP  (AB)  ,  686  for  RKKI' 
a3C).  and  189  lor  PnKP  (AB)  . 

Typical  example.s  of  histograms  of  summary  residual  data  are  shown  in 
Fig.  711.2.1,  III  tbi.s  figure  we  have  also  plotted  (he  data  from  deep 
evciit.s.  A  compari.son  suggests  tliat  reading  errors  are  significant 
f.sp.  r  iaiLv  ior  ai  i'ivals  from  shallow  events;  note  that  the  luimbei  of 


late  readings  is  reduced  in  the  data  from  deep  events.  For  the  core 
phases  and  for  P  at  distances  larger  than  85°  (P2) ,  the  early  parts  of 
the  histograms  for  all  data  and  for  the  data  from  deep  events  overlay 
quite  well.  This  means  that  the  upper  mantle  model  of  PREM  is  consis¬ 
tent  with  the  ISC  depth  estimates.  The  P  data  at  distances  smaller  than 
85°  (PI)  are  different  in  that  there  are  anomalously  many  early 
arrivals  from  shallo<j  «.v/ents;  thi.s  may  also  explain  the  relatively 
large  variance  of  these  data.  It  is  possible  that  one  begins  to  see 
here  the  effect  of  subduction  zones,  since  many  of  the  events  occur 
within  these  zones. 

Fig.  VII. 2.1  also  shows  that  there  is  a  significanr  mean  residu.il  Icf: 
in  the  data.  This  is  especially  clear  for  the  core  phases,  and  we  can 
infer  their  relation.  If  the  sampling  by  summary  r:iys  is  reasor;ably 
uniform  and  if  nonlinear  effects  can  b(^  neglected,  thtui  for  any 
particular  phase  the  mean  residual  represent';  the  ifiects  ol  dil- 
ferences  between  PREM  and  the  spherically  averaged  earth,  and ''or 
systematic  reading  errors.  The  PeP,  PKP  and  PnKP  mean  residuals  for 
summary  rays  in  the  same  r-^y  parameter  interval  are  expected  to  follev.- 
a  linear  trend: 

^(PnKP)  ■=  +  niTc.  n  1,2 _ 

where  ST^p  =  5T(PcP),  and  5T^  represents  the  residual  aftei'  one  pnssa.i 
of  the  wave  through  the  core;  both  the  velocity  stricture  and  the  core- 
mantle  boundary  level  may  contribute  to  the  residual.  A  relation  of  •  In¬ 
form  (1)  can  be  discerned  for  the  phases  with  ray  paraiiietf  rs  above 
'4  s/d,  but  surprisingly,  PeP,  PKP  and  PKKP  in  the  my  parameter  range 
2-3  s/d  do  not  follow  a  linear  trend.  One  possible  explanation,  now 
under  investigation.  Is  based  on  the  fact  tliat  PeP  -it  small  distances 
is  weak,  and  known  to  be  often  unobservable.  It  is  therefore  po.ssible 
that  PeP  (and  possibly  PKKP)  is  observed  primarily  in  c i rcums tance;  of 
relatively  strong  focusing,  with  an  accompanying  phase  dti.iy. 


It  ir:  also  of  interest  to  iioti'  that  tlie  variance  of  the  PKKP  data  is 
not  much  larger  thaii  that  of  PKI'  in  the  .same  ray  parameter  range  (?-3 
s/d}.  It  is  convenient  to  plot  the  variance  of  the  various  phases  as  a 
funr  tion  of  tht  ir  seii.siti vi  ty  to  varintioTis  of  deep  earth  sti  ucture. 
Hero  we  give  such  a  relation  be  tween  the  variance  of  the  data  c-p  and 
t'ne  variance  ot  core -mant  1  f-  hound.ary  t  opograpiiy  For  PeP: 
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and  for  PnKP  if  the  piu' tvu't>  i  ons  6r  in  tlie  sampling  poirt.s  r^  are 
uncr  rrel a  ted : 
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Here  y\  =  r/v,  and  a  .superscript  refers  to  the  top/ho  t  ’  ()nis  ide  of  thf 

i.)ou’.idarv .  In  Fig.  VII  2.3  the  variance  of  tiie  data  subsets  is  plotted 
following  equation  i 2 ) .  One  inference  from  this  figure  is  that  the  PKFl 
dat.i  Lmpl\'  a  re  I.  a  tivc  i  smooth  core-mantle  l)oundarv  on  ,i  large  .scale; 
for  i  ■  1 '  rs  •  ra  t  i  V,  purpr.je.s  tin;  (■\[u-cted  travel  time  vari.ince  for  - 

i  is  siiown  !  Fig  V I  1  . .  3  .  Anotiier  inference  i  ,s  tl'u.t  models  ol 

l.irg.e  -  sea  1 C'  later.il  '.’ariation  of  dc-ep  earth  structure  can  tspl.iin  onl' 
a  re’.arivclv  small  p.-.r;,  of  the  d.il  a  variance. 
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Fig.  VII.2.1.  Histograms  of  travel  time  lesiduaLs  of  'summary  rav' 
data.  Ray  parameter  ranges  are  2-3  s/d  for  the  BC  branches  of  PKP  and 
PKKP  and  for  PcPl,  3.5-4  s/d  for  PcP2  and  >4  s/d  for  tho  AB  brauclies  ot 
PKP  and  PnKP  and  for  PcP3.  Distance  Interval  for  PI  is  65-85°.  for  P.' 

85-95°.  :  all  data;  .  :  data  from  deep  events  (  >'i()0  km 

depth) . 
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Lk&j — V 1 1 . 2 . 2  ^  Mean  'summary  ray'  Lravel 
PnKP,  n  >  ?.  o:  2  <  p  <  J  f-./d;  •:  p  >  /, 
linear  fit  to  t:h('  data  v;i  th  p  >  /(  s/d. 


t  ime  residuals  of  I’d’  and 
s/d.  The  dotted  Line  is  a 
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of  \’('r  t  liSoa 


Thi::  is  the  thii'd  and  final  report  on  iiiodolling  of  Lg  wave  propagation 
in  tho  Central  Grahen  of  the  North  Sea  in  an  attempt  to  explain  the 
ver;;.  strong,  atienuation  of  the  l.g  wavetrain  obscrv'ed  in  this  area.  In 
the  .first  report  (Maiipin,  1987),  we  presented  the  modelling  method  and 
son;r  preliminarv  te.st:.  The  bulk  of  the  modelling  re.sult.s,  i.e.,  the 
reflection,  and  transji!  i  .ssi  on  matrices  for  Rayleigh  and  Love  type  Lg- 
modcs  propagating  at  ■.  right  angle  or  at  an  oblique  angle  across  a 
graben  model,  were  presented  in  the  second  report  (Maupin,  1988).  A 
first  interpretation  of  the  matrices  showed  that  on  the  average  over 
manv  Lg  wavv-rra  in.s ,  80t  of  the  incoming  Lg  energy  remains  in  the:  Lg 
v.nr.'c  after  prop  igation.  ac  re  s  tlie  graben  model. 

The  transmis.s  ion  affects  d  i  I  {<•  ron  1 1  v  incoming  Lg  wavetrains  wi  tl 

different  modal  contents.  Th(  rel.ilive  amplitudes  of  the  different 
modes,  which  dencnul  on  tlioir  e.witation  by  the  .seismic  source,  tis  weli 
as  tiieir  phase  (,ii  f  leronces  ■..lien  reaching  the  Graben,  which  vary  with 
epicontrai  dista.nce,  define  the  modal  concent.  In  order  to  eMploit;  mori. 
coiiipletciv  the  tivtr  smi  .ss  Loii  matrix,  we  analvze  bore  its  effect,  on  Lg 
v.’a'.-f"  raln.s  from  dilfcrent  sources  at  differc-nt  d  1  .s  t..r.ec-is  from  the 
Grai.'tn,  Since  the  r  ransm  1  s.s  i  .>n  of  R.ivleigh  and  i.ove  waves  at  a  right 
angle  or  ac  an  ('bliqui  angle  .icros.s  the  model  have  been  found  very 
similar  in  the  secoml  repor;  ,  tje  concentrate  fur  analvsi.s  here  to 
Rayleigh  waves  propagating  at  right,  angles  across  the  structure.  On  :  n. 
other  hand,  the  t  ransi'i  i  s.s  i  on  acro.s.s  three  variants  or  the  Central 
Graben  model  used  in  tiie  [u'f 'cious  i'  |>orts  (and  now  called  model  1)  arc- 
al.so  examined,  to  accoimt  for  po>.-.ilile  I)  1  ock  -  fau !  t  ing  of  the  Graben 
mr.rgln  (models  2  and  Fig.  VII,  ’.11  or  roughnes.s  of  the  sc  diment - 
baseiiiCnt  interface  (model  4,  Fig,.  VII.  1,1  i.  We  at.so  extimine  the  pha.si 
stability  with  i.ci'iod  .'f  the  t  r.iasiiu  l  t  ed  wavetrain. 

Aftei  inspect] or  of  the-  results  lor  different  souices ,  we  retain  three 
typii-al  cases  for  discussion:  an  explosLon,  a  strike-slip  cvirClujuake 
with  a  fault  trace  a'  / from  the  syii'm*  t  rv  dinet  ion  of  tlie  Graben. 


and  an  earthquake  which  occurred  on  the  western  flank  oL  the  Viking 
kraben  on  29  July  1982  (strike;  100°,  dip;  63°,  slip;  -170°,  after 
Havskov  and  Bungum,  1987),  for  which  we  study  Che  waves  travelling  due 
east  perpendicularly  across  the  Viking  Graben  and  to  Norway.  We  use 
explosions  with  4  different  focal  depths,  ranging  from  0  to  3  km,  and 
earthquakes  at  7  different  focal  depths  sampling  the  whole  crust.  TIil; 
distances  of  the  events  from  tlie  Graben  are  taken  rangin,  fro;r.  G  to 
1000  km,  with  a  step  of  10  km,  providing  a  good  sampling  of  possiblt 
phase  shifts  between  the  different  modes  when  reaching  tiu‘  Graten. 

These  events  do  not  intend  to  model  a  complete  or  realistic  situation, 
but  to  provide  an  oversight  of  the  effect  of  the  Graben  on  tliff.ront  I... 
wavetrains.  We  recall  that  Gregersen  (1984)  used  many  ta  rtliquakc  s  ii; 
his  study  of  the  attenuation  across  the  North  Sea  Central  Graben .  ar..; 
pointed  out  that  the  effect  does  not  depend  on  the  source. 

The  total  energy  transmission 

For  each  source  type,  depth  and  distanct*.  wc  c.il  cu  i,i ;  e  the  aii'our.t  of 
total  energy  contained  in  the  Lg  wavetr.iin  before  ami  aitir  propetgat  ion 
across  the  Graben.  This  total  energy  Includes  tiie  eiu  rgv  contained  iii 
the  whole  crust  for  the  11  Lg  modes.  For  each  .sourr<  nie')i.,n  i  .sm  .,nd 
depth,  the  results  are  summari:;c;d  in  a  lustogram  of  tin-  •  .  ■  s  i,,;. 

ratios,  which  illustrates  how  tlieir  values  vary  for  diftciviit.  source- 
graben  distances. 

One  of  these  histograms  is  shown  on  Fig.  VII.  2, 2.  1'  dispiay.s 

values  of  the  transmission  ratios  across  model  1  foi-  ait  :  ,  '.m'.  i,  ;  ra  it: 

excited  by  a  Viking  Graben  eartliquake  at  13  km  local  dcprli.  Th^ 

distribution  is  well  peaked  around  transmission  ratios  o(  >:n  . 
lli.stograms  for  otlu-r  focal  depths,  source  inechanisii"  oi  i  ■  K  i ,i  (  'a  i  ■. 
similar  in  shape,  with  a  filight  shift  of  -10"  ('u  a  .pi  a.inia,  ..l-i.s,-  to 

the  surface.  Fig.  VII. 3. 3,  wheta,-  incident  .nid  mean  :  i  nn-S"' i  ;  ;  i-d  onrr('i< 
are  plotted  as  a  function  of  source  type  and  depth,  also  tistifics  tiiat 
in  the  large  majority  of  caacn.  80%  of  the  incident  ener,  ■■  is  trans¬ 
mitted  as  an  Lg  wave  across  our  models  of  the  Centr.il  c;ra!a  i,.  Tl.. 
remaining  20%  is  converted  to  Sn  or  other  S  waves  propagating  in  the 
mant le . 


Thif-  resvilt  is  lu  agrooment  with  t;he  fincliiigs  in  report  no.  2,  and 
shov.:5  that  the  total  energy  transmission  ratio  is  only  slightly 
dcpi:.dont  on  the  soun  e  mechanism  which  has  excited  the  Lg  wavetrain. 

The  surface  energy  transmission 

Since  .surfiice  waves  liave  their  energy  distributed  with  depth  dif¬ 
ferently  from  one  mode  to  the  other,  their  total  energy  does  not 
directly/  indicate  iiow  much  of  the  energy  is  confined  ciose  to  the 
surface,  or  equivalently  the  surface  displacement.  More  in  agreement 
with  what  can  nt^tuall;.’  be  measured,  we  therefore  also  analyze  the 
ratios  of  tran.sii'i  t  t  od  over  incident  surface  energy.  This  surface  energy 
is  the  energy  of  the  whole  hg  wavetrain  measureci  on  the  vertical 
coii.ponont  of  o  sc  i ,  xho  ratios  of  transmitted  ovei'  incident 
maximum  vertica’  isp  1  .icemeiit.  at  the  surface  were  also  calculated,  bnf 
are  not  ri.,scus.svd  her.  .since  the  more  global  rhnract,;r  of  the  energy 
make."  it  prio.-j  a  iin' re  st..hle  .'luantity  for  estimating  the  attenua- 
rioit.  Wo  do,  however,  observe  a  high  degiwe  of  similiirity  between  the 
ratio.s  in  energy  and  in  maximum,  displacement. 

The  pattern  of  surface  energy  transmission  is  verv  different  from  the 
pattern  of  total  energy  transmission,  'liireo  typical  histograms  of 
surf.ice  energy  r  ransm  i  as  ion  iwr  i.'s  for  difiin'cnt-  source  -  graben 
distances  are  disjilayed  on  Fig.  \  1  I  .  3  . 'i  .  and  incident  and  mean 
!  r  ciMii  i  t  tc  d  surface  .'nergies  a.‘;  r.  fuin'tioii  of  .source  et.'pth  are  plotted 
in  i : g .  VII. 3. 5  and  VII. 3. 6  for  different  sources.  In  order  to  indicate 
the  dispersion  in  the  t  ran.smi.ss  i  on  due  to  rlianges  in  tlie  irodels,  w. 
plot  the  maximum  and  minimum  v-aluc.';  oi  the  mean  transmitted  energies 
calculated  with  aiodels  1  to  1  .  Mo  ruK-  applies  as  to  wiiich  model 
usually  gives  the  lower  or  higher  v.alue. 

Fig.  VII. 3. 4a  is  a  typical  histogram  for  explosions  c.r  very  shallow 
ear '  nquakes ,  for  whic  h  the  surface  transmi  ssion  ratios  arc^  smaller  than 
tile  total  energy  tfansmis.sion  ratios.  For  these-  sourcc.s,  a  large  part 
of  the  tot.il  T.g  energy  is  confined  close  to  tlie  surface  of  the  model, 
malu.v  in  the  sedimentary  layer,  liefore  reaching  the  (iraben.  Crossing, 


the  Graben  shifts  part  of  the  energy  deeper  in  tlie  c mist  by  redistri¬ 
buting  the  energy  more  evenly  among  the  different  Lg  modes.  Thi.s  effc-ct 
amplifies  at  the  surface  the  global  loss  of  Lg  energy.  The  surface 
signature,  of  the  Lg  wave  i.s  therefore  decreased  by  a  factor  of  0.6  Zj 
0.7S  in  terms  of  mean  amplitude  of  the  .signal. 

Some  earthquakes  with  certain  focal  mechanism.^  or  Located  at  tl.e  botioin 
of  the  crust  excite  more  ev'enly  the  different  modes  oi  tlie  Lg  wave.s. 
This  is  the  case  for  mid-crustal  or  deep  strike-slip  ear thqu.akes ,  for 
example,  which  have  rather  well -peaked  transmission  ratio  distr ibution.s 
(Fig.  VII. 3. 4b),  similar  to  those  for  the  total  energy,  and  mean  value.s 
of  surface  energy  transmission  around  80%  (Fig.  VTl  in  tha“  case, 

the  mean  surface  displacement  is  decreased  by  a  factor  ol  0.9  after 
cro.ssing  the  Graben,  and  this  directly  accounts  foi  tliC;  total  loss  o' 
C'nergy  in  the  Lg  wa-.'etrain. 

Other  types  of  mid-crustal  earthquakes,  like*  the  mid- crustal  Vif.ing 
Graben  earthquake,  excite  primarily  the  Lg  mode:;  having  i!ieir  energy 
confined  in  the  middle  of  Che  cru.st.  The  surface  energy  before  reacliing, 
the  Graben  is  thus  small  compared  to  the  total  energy  involved  in  th* 

Lg  wavetrain.  By  crossing  the  Graben,  the  energv  is  redistributed  anc>ng 
the  modes,  and  some  energy  is  thereby  shifted  from  tlie  mi  Idle  of  tl.' 
crust  towards  the  sedimentary  layer  and  tlie  surface .  The  net  et  ioct  :  ,s 
an  increase  in  surface  energy  (Figs.  VII. 3. 4c  and  V!1.3.')',  despite  ti." 
decrease  of  total  energy.  In  that  case,  an  ita'rease  of  1.  .  can  iic 
expected  for  the  mean  amplitude  of  the  recorded  Lg  t  r.i i  n . 

The  total  energy  transmi.ssion  ratios  have  shown  th.it  i  1k  doho  ■  a 

rather  energy-proof  barrier  (only  20%  of  the  energy  lenk.s  into  the 
mantle).  On  the  other  hand,  the  surface  energv  ratio.s  .slur.-  that  the 
crustal  thinning  of  the  Central  Graben  causes  inipori.aul  transfer.s  of 
energy  among  the  different  vinits  of  the  crust.  Comparing  the  .suiface 
energy  curves  (Figs.  VII. 3. 5  and  VII. 3. 6)  with  the  total  energv  ones 
(Fig.  VII. 3. 3).  we  see  that  the  surface  energy  curves  are  very 
different  from  the  total  energy  ones  before  propagation  across  the 
C.’raben  (filled  symbols),  but.  much  more  similar  aftcrwr.rd;.  I'c.pen 


symbols).  The  L.'raben  has  redistributed  more  evenly  v\;ithin  the  crust  tiie 
total  energy  in'.’oL^’cd  in  tlie  Lg  wavetrain,  and  the  surface  energy 
reflects  better  the  total  amount  of  energy  contained  in  the  whole 
crus  t . 

Propagation  across  our  Graben  models  leads  to  some  Lg  amplitude 
variations  at  tne  surface,  though  limited  in  size  and  both  positive  and 
negative.  They  are  very  different  from  the  factor  0.25  to  0.5  actually 
observed  in  the  North  Sea  Grahen  area.  Moreover,  the  preferred  focal 
depth  of  seismic  events  in  the  North  Sea  is  very  often  around  15  km 
(liavskov  and  Ruiigum,  l'!87),  whicti  would  bia.s  our  transmission  ratios 
towaids  their  higtie.s;  ■.•alne-;.  Our  modelling  wc.uld  at  the  most  explain  a 
factor  of  2  between  the  attenuat.  ion  of  Lg  waves  produced  by  explosion.^ 
and  Lg  wa-'es  preJured  by  earth<.|n.ike.s ,  hut  can  in  no  case  explain  the 
general  and  strong  a !  •amuat  i on  observed  in  this  area. 

Coherency  of  the  phase  with  period 

The  previou.s  ca cul  a  t  i  nns  have  been  made  at;  a  single  frequency.  The 
phasc'  bt.'havior  of  the  waves  as  a  function  of  period  is  a  key  cilement  to 
the  effective  bulld-ui->  of  a  wavetrain.  If  rapid  variaticins  are 
observed.,  interfcri  Ufa  betveen  noighborinr,  ptriods  might  destroy  tl';e 
t  .'‘ai In  ord<M'  m  cluu  I,  the  stabilit.y  of  t.ho  pha.so  a.o  a  function 
o'"  p  -riod.  we  n.ov;  co:;';  ire  t  -ie  ph  i.-a  .s  of  thi-  hr  wave  mode.s  propagat  itif, 
out  c' I  the  Geiitrol  c'r, at  ?  ut  i  gjiivu- i  nr.  period;;,  1,0  and  1.02  s 

We  calculate  the  t  imn.smi  s.s  i  on  matrices  at  the  periods.  For  the  same 
series  of  sources  and  souree  -  grabeti  distances  a;;  e.-.rlior  in  this 
report,  we  calculate  tlie  phase  of  eacli  l.g,  mode  propagating  out  of  the 
Grahen  at  the  2  perioci.s.  We  must  note  that  a  mode-  prc'paga  t  i  ng  out  of 
tile  Graben  or  i  g,  i  n.i  f  r  .s  from  the  combin.it  ion  in  the  Ginaben  of  different 
modc'S  initially  vcscitcd  by  the  source.  Its  phase  thus  depends  in  a 
complicated  way  ,ni  tlu  phase.s  of  1  liese  modes  wlien  they  enter  the 
structure.  We  subtract  from  the  lot.al  phase  the  pure  propagation  pliase . 
i.e.,  the  integral  over  horizontal  ilisianct  of  the  mode  local  phase 
slowiKss.  By  u.s  i  ng,  a  pliase  free  of  pure,  prop.ag.i  t  i  on  effect,  the  phasf 
diff c  rc  nee  t)otwef;i  r  hi  mode.s  .it  t  l.c  2  different  periods  is  actually 


measiired  at  the  arrival  time  of  the  modi;  predicted  by  its  group 
velocity . 

On  Fig.  VII. 3. 7  are  displayed  three  histograms  of  the  phase  dif¬ 
ferences  for  the  different  modes  and  different  source-graben  distance5; 
after  propagation  across  model  1.  Due  to  the  unknown  but  certainly  poor- 
accuracy  of  the  transmission  matrix  phase,  which  is  influenced  by  the 
zoning  of  the  model,  we  cannot  use  these  histograms  verv  quanti¬ 
tatively.  Even  if  in  the  second  one  large  phase  differences  occur 
rather  frequently,  cases  a)  and  c)  testify  that  the  phases  of  the  modes 
are  not  systematically  random  after  crossing  the  Graben,  and  therefore 
cannot  give  rise  to  a  generally  strong  attenuation  of  thc'  wavetrain  by 
destructive  interference. 

Conclusion 

The  investigations  presented  in  this  report  cunfirin  the  cunlusions 
already  drawn  in  the  second  report. 

Our  numerical  modelling  of  Lg  wave  propagation  in  a  simj'l’fied  model  ot 
the  North  Sea  Central  Graben  doe.s  not  predict  the  i-u'.a  re  attenuation  uf 
the  wavetrain  actually  observt>d  in  this  region.  On  :  lie  contrary,  the  l.i, 
wavetrain  appears  very  robust  when  cros.sing  a  .tone  where  its  waveguid- 
is  strongly  deformed. 

Since  the  large-scale  geometry  of  the  Graben  tails  to  oxilain  tiie 
observed  data,  we  suggest  that  future  work  explore  al:-ernitive 
explanations  for  the  observed  attenuation.  .Scattering  l)y  .  Ij  or  ii) 
basaltic  intrusions  in  the  lower  crust,  extensive  fault  ii."  issocL.i'.c 
with  intra-fault  weak  material,  or  more  rheoloj-ieal  aspoe  .s  niiy'it  he 
goot!  candidates. 


V.  Maupin,  i'os  td<-c  t  o  ra  ;  . 
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Fj  F.  ■  V'l  I  ■  3  ■  I  ■  Model.s  of  the  North  .Seo  (.'cMttr.t!  G'raiu'ii. 

Model  1:  Full  line  model,  used  in  the  previous  reports 

Models  2  and  3:  Block- faulted  models 

Model  4:  The  same  as  Model  1  with  perturbations  of  the  sediir.en,  - 

basement  interface  represented  by  a  dotted  line. 
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F’ig  ■  VII ,  3  ■  2  ■  H  i  fi  f:o);r;ira  f)f  ciie  eiu'rgy  Lransmiss  inn  ratios  for  n  Vik 
Graben-type  event,  at  1  a  kiti  fo<-al  depllt  anti  di.sLance.s  from  '  iu-  (•rabcii 
ranging  from  0  to  ]  (K)(i  km. 


energy 


Fig.  VII  .3.3.  Total  energy  in  the  Lg  wave  betore  (filled  .symbol)  and 
after  propagation  across  model  1  (open  symbol),  as  a  fvmctlon  of 
source  type  and  depth.  The  energy  scale  is  only  relative  .since  no 
physical  source  size  is  included  in  the  modelling. 
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Transmission  ratio 


Fi  p  .  VII  ■  3  ■  4  .  Histogr.-iin  oC  surface  energy  transmission  ratios  for: 
a)  an  explosion  at  ihe  surface;  b)  a  strike-slip  event  at  IS  km  focal 
depth;  and  c)  a  Viking  Graben-type  event  at  15  km  focal  depth  and.  for 
all  rases,  di, stances  from  the  Graben  ranging  from  0  lo  1000  km. 
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P’ig.  Vll.3.5.  Surface  energy  in  the  Lg  wave  measured  on  'he  vert.ici] 
component  before  (filled  symbol)  and  after  propagation  acrof;s  Lie 
Central  Graben  (open  symbol),  as  a  function  of  source  type  and  liepth. 
The  minimum  and  maximum  values  of  transmitted  surface  encigics  ;  verag  d 
over  aifterent  source - graben  distances  for  the  four  Grabei'  models  art 
represented.  The  energy  scale  is  only  relative  since  no  physical  source 
size  is  included  in  the  modelling. 
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Phase  diMerence 


!'  i  i' .  V  i  I  .  i  .  7  .  H  L.s  t  oc.ram  of  the  [)(iase  difference  bf-[w<-cn  1.,;  iiiodis  ,1 
pcriod.s  ]  .  fl  .ind  1  .0'/  s  aftt-r  ['I'opnp;;!  t  i  on  ncro.s.s  model  1  .  lor:  .1)  ni 
I  v[);o.sinn  .'It  the  .surface;  b)  a  Vikinp.  (.’raben- 1  ype  <  vent  luii  f 

(le()iii;  and  c)  a  Viking  Craben-type  event,  at  ]  !)  km  local  d.  [dli,  .asl, 
a  1  ca‘.e.s  ,  distinces  f  10m  tlie  graben  ranging  from  d  to  ll'v'  km. 
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:* _ The  August  8.  I9PR.  More  Basin  eavthnuake:  Observed  Kround 

raotiiM.s  and  inferred  so'irce  parameters 

An  i-artrli(;  ;:ii;!.'  (■>!'  niaj’.!) •' tudo  >irouiiil  b.2  occurred  in  the  Motf  H.isin  o’\ 
Auju'it  8,  ]d88,  with  tremors  lelt  ovei"  most  ot  soutiiein  and  cc'ntral 

Nor'.'uv.  The  eai  thquake  v;as  tlie  l.u'gest  one  in  tiie  rej’.ioi)  tor  at  least 
iO  '.ears,  vitii  a  focal  mcchanisin  solution  that,  indicates  thrust 
faultini-;  aloni;  a  KNE-SSW  .strikinj;  fault  plane,  i  it  response  to  E-iv 
compress  i  onal  stress.  The  seismic  moment  was  of  tlio  order  of  10^^  Nit 
wi  til  indications  of  a  scaling  consi.stent  with  an  lo-square  source  riodc  I 
A  major  .■-.ource  i  !'  nneertaint  v  in  this  analysi.s  is  tied  to  anela.stic 
attt  caati  an . 

Background  seismic ity 

The  seismicity  of  tills  part  of  Norway  and  rlie  N'orvtegian  Continental 
Shelf  is  shown  in  Fig.  VII. 4.1,  wlio.re  three  diffi'i'ent  time  periods  hre. 
been  plotted,  with  different  symbols.  The  map  indicates  a  reasonably 
gco-.l  eorrelatioi;  he:,vt<,fn  s.'i.smicity  and  rcgioijal  geological  features 
-Mch  as  fault.s,  fault  c.ovu  .s .  fracture-  zones  and  p,rabon.s ,  and  thf-re  .iia 
also  indication.'.;  c'f  the  seismicitv  following  the  continental  margin. 
The  m.'ip  moreover  show.s  a  certain  hia.s  between  the  different  t  i  me 
periods  in  that  the  areas  .south  ot  C.l''N  obviously  arc  btitrer  co”.:red 
tf.-rm.;  of  n:  i  crose  i ';;ri  i  c  .survi:  i  !  I  aiu  >  dui  inj',  the  ! '-tSOs  (Buiigum,  1 )  . 

Epicenter  location 

The  location  of  'he  August.  E,  1988,  e.irthquake  in  thi.  More  l-u.ii.  i  .s 
shown  in  Fig.  i  :  .  d  .  1  .it  hi  /®N,  .slightlv  west  of  most  e.mlt.  : 

events  in  this  area,  hut  still  e.ast  ot  the  prominent  F.trnc- -  She- 1 1  and. 
Escarpment.  The  ear'i.  heuake  was  widely  recordi-ct  on  sv  isniic  i  nstruiiH'ni  .- 
throughout  northern  Europe  and  tlie  entire  world.  Nearly  80  seismic 
phases  have  been  reported,  all  within  the  distance  range  of  100  to 
I  lot)  km.  F.xper  i  nw'iii with  loc.-iLing,  t  lie  event,  with  suhs'eis  of  those  d.,: 
insured  the  con-,  i  ti-ns  and  reli,d>ility  of  tlu-  location.  No  n-Liihl'' 
depth  estimate  i  vc-l  ,<vailah]e,  however  fHan.son  et  ,i  1  . 


Felt  effects  and  magnitude 

This  earthquake  was  widely  felt  throughout  much  of  central  and  south. in 
Norway,  along  the  coast  from  Stavanger  to  Mo  L  Ran.i .  as  veil  as  in 
southeastern  Norway  and  in  Sweden,  Responses  to  quetitionnaires  sent  out 
ny  the  Seismological  Observatory  in  Bergen  give  felt  radii  of  about  100 
and  440  km  for  intensities  IV  and  III,  respectively.  In  rising  rtlation- 
ships  developed  recently  between  felt  area  and  surface  wave  magnitude 
Mg  (Muir  Wood  and  Woo,  1987),  this  results  in  Mg  values  of  3.2  and  5.3. 
respectively.  In  comparison,  an  Mg  value  of  5.1  +  0.24  has  been 
computed  by  N.N.  Ambraseys  for  this  earthquake,  while  NORESS  data  ba-.'e 
given  an  value  of  5.2.  This  magnitude  makes  this  earthquake  the 
largest  one  in  the  rc'gion  for  at  least  30  years,  possibly  even  the 
largest  one  since  1895  (Hansen  t't  al,  1988;  Bungum  and  Seines.  1983; 

Focal  mechanism 

The  .sense  of  faulting  for  this  C'.arthqu.'ike  w.ns  explored  througli  he  ■ 
of  the  direction  of  vertical  motion  of  about  30  of  the  iirst  arriviin; 
P-pliases  for  all  available  recordings.  A  focal  mechanism  solution, 
using  this  approach,  is  givtui  in  Fig.  VII. 4. 2.  where  a  combination  of 
local  and  teleseismic  data  helps  in  constraining  the  nodal  plane.s.  Fvoiii 
the  graph,  the  faulting  parameters  for  the  two  planes  are  strike  20°E, 
dip  46®  and  rake  (slip)  116®,  and  strike  165®B,  dip  50®  and  rake  66®, 
respectively.  The  solution  leaves  an  ambiguity  as  to  which  of  the  cwi 
planes  is  the  faulting  plane,  hut  in  either  case  this  solution  gives, 
reverse  mechanism.  From  the  geologic  data,  however  i see  Fig,  Vli  9.1 ; 
we  find  that  the  preferred  fault  plane  for  this  earthquake  i  .s  th.c  oik 
striking  20®E.  It  can  be  seen  from  Fig.  VII. 4. 2  that  this  nor  tin  as  t  e  r !  .■ 
striking  plane  is  only  constrained  by  the  stations  whose-  azimuths  var 
from  about  40®  through  about  90®,  These  are  the  stations  nt  the  SEISN  1< 
network  and  ARGESS.  The  sense  of  first  motion  changed  from  dilatation 
to  compression  through  the  middle  of  this  netv.’ork  of  .stations,  tillowii  a 
tor  a  wel 1  -  cons t rained  fault  plane  solution  (Hansen  et  al.  1988), 

Observed  ground  motions 

An  earthquake  of  magnitude  5.2  naturally  causes  most  convent  ion. i ! 
seismometers  within  regional  distance  ranges  to  saf'u'uic  iTiclip'n  1 


record inr,:.  have,  howe'.'er,  in  tlic  present  case  betn  obtained  at  throe 
sites;  (li  Moldo  (MOi.j  accelerometer  site  withiji  the  SEISN’OR  (Northev! 
Norw.iy)  network  (288  k.m)  ;  (?)  5>ulen  (SUE)  acce  I  orometer  site  within  the 
Western  i,’orv/ay  network  (818  km);  and  (3)  NORESS  HE  (high  frequency)  and 
IP  (intermediate  period)  eltinent.s  (578  km).  Since  all  of  these  moreovor 
yield  broadband  recordings,  they  become  especially  valuable  in  terms  of 
infe ’"encos  about  source  paramei.e  rs  (NORSAR  and  Risk  Engineering,  I  uc  .  , 
19881 . 

Observed  source  disp 1 ..  ement  spectra  for  these  stations  are  shown  in 
Fig.  VII.  i*.  3.  where  t''.e  time  series  were  rotated  to  yield  the  radial 
(R)  and  tiie  transver.S'  (T)  components.  The  data  are,  corrected  for 
svstem  re.spor..s.'  .  including  .special  proce.ssi  r;g  of  the  accelerometer 
data!,  and  converted  from  a.celeration  to  displacement  for  Molde  and 
Su;e:i  at.d  from  \'c-loci“''  to  di.spl  acement  for  NORESS  (including  a  careful 
bandpass  ’.  ilteriug  in  both  ca.se.s).  Energy  spect.ra  are  then  estimatf  d  as 
a  basis  for  the  plotted  tii.sp  1  .nceii’cnts ,  with  a  time  window  covering  29 
seect id..s  of  the  v.iv. 

'..'hat  i  .s  .seen  f  la  :ii  Fig.,  Vn..,i  if  that  t  lie  ob.si'-r'.s-^d  d  i  sp  1  .areiiient  s  fal; 
(jff  vith  frequenc'.'  .it,  ■.  i\d  i,-  not  ■•’.■rv  diffeu'ent  froiii  (as  indi('ateci 
!iy  straight  liiii.-)  .At  i-.iidcr  f  lu  queue  ies .  the  .slop"  -bereases 
some'what,  possiblv  ini'uer.cc-d  by  noise,  and  at  low  ■.  a.'t  t'nc  i  es  it 

should  be  kept  in  mind  that  the  .spectra  arc  certainly  affected  by 
noise.  The  filtci.s  usa-d  in  proce.'U'.ing  these  data  have  been  cUfined  at 
lower  cutoffs  at  0.20  'Az.  for  Molde  (where  quantization  noise  also  niav 
tiav"  been  ,a  problem)  .cid  S-il«-n,  at  0.1''  Mz  feu'  NORESS  IP,  ,ind  at 
0.80  Hz  for  NORESS  HE, 

Corrected  ground  motions 

In  order  to  be  ahli  to  i-oiiqiari'  the  ob.sirved  g.r<>und  motion  d  i  spl  aceii'c  n ; 
■Spec  t  I'.i  more  coir.'f-ni  c  lu  ]  y  ,  v.e  liave  i-or  rec  t  id  all  of  them  lor  the 
effe'  ‘  s  of  geoiiu  t  rival  .spreading  .iiicl  aiielast  ic  .itlenuation  hack  to  .i 
reference-  distar-.  of  10  km  1  rom  (lie  source,  with  results  a.s  shown  in 


The  correction  used  for  geometrical  spreading  has  been  the  commonly 
used  model  by  Herrmann  and  Kijko  (1983)  in  which  tliere  is  a  change  from 
spherical  to  cylindrical  spreading  at  a  distance  of  100  km; 

R  <  100  km 

G(R)  = 

0.01(R/100) R  >  100  km 

For  anelastic  attenuation,  we  have  for  test  and  sensitivity  purposes 
used  two  very  different  models,  by  Kvanime  and  Havskov  (1938): 

Q  =  120  •  fl  ■  1-  (Model  1  ) 

attd  by  Sereno  ft  al  n988): 

Q  =  560  ■  fO-26  (Model  2) 

The  first  of  those  has  been  developed  from  sptci  rtil  r.it  h'  .n  d  c,  da 
decay  methods  based  on  data  typically  within  a  100  100  kit  distance 
range,  while  the  second  has  been  developed  from  a  s i mu  1 1 .a  o'.'S 
inversion  for  seismic  moment  and  apparent  attenuation  l;)ast.  <i  on  uat.i 
between  200  and  1400  km.  The  first  model  covers  2  - 1 'i  H.; ,  and  the  set.''!' 
1-7  Hz. 

Kven  with  these  differences  in  mind,  it  is  no:  ol>viou'.  wiii<-h  oi.e  of  'hi 
two  models  will  be  most  appro)>r iate  in  the  pr*  .S'-nt  '  i  i.li  .‘■■.rl, 

large  differcnce.s  in  Crtiquency  .'.oii.slt  ivi  t  y  o1  t  l.t  !.'<  th  is.  !  io  '.'e  r  ,  i  ; 
is  understandable  that  the  effect.s  of  ;  he  parli  e  or  ri.  ct  i  on  i  will  bo 
different,  as  shown  by  Fig.  VI  (  .  i  .  4  It  t.s  useful  i,.  re  to  note  ;!;at;  ’i 
the  observed  and  corrected  spectra  have  .slopes  proportional  to  and 

1^,  respectively,  then  the  following  frequency  sen.s  i  t  i  i  ; of  the  0 
model  will  be  required  (Q  »  Qgf9): 

vQ„(i -7)lnl0 


n 


1  -  log[l  + 


1 


whera  v  Ls  wavt  '^eloi  ity  and  R  is  distance.  This  relation  sliows  that  if 
one  rer|iur,..s  th.  ohs,  i  ved  .slope  to  be  maintained  after  corection  <  S 
>),  then  Q  must  be  directly  proiiort  ional  to  frequency  (rj  --  1). 

.'•’roll.  Fig.  VII. .'  i;  i  .s  seen  that  the  ob.served  spectral  slope.s  are 
reasonably  close  to  for  most  of  the  data,  at  least  in  the  l-o  Hz 
range.  This  slo-pe  i  •■i  therefore  more  or  less  maintained  through  the  path 
correction  when  using  the  Model  1  attenuation  (»j  =  1.1)  as  shown  in 
Fig.  VII. -1.4,  '.•.hich  in  turn  gives  indications  of  a  source  model  close 
to  t  hie  standard  Bi'une  model  (see  also  Chael  and  Kromer,  1988).  The 
Mode- 1  2  attenu.r  ion,  -.n  the  othei  side,  is  more  difficult  to  reconcile 
•v’ith  this  parti  .'ul  ,1 set  ol  d.»ta. 

Source  displacement  spectra  and  seismic  moment 

From  Che  observ.  i!  di ';pl  aceinent  sp<  ctfi  in  Fig.  Vll.4.3  (or  from  tiie 
corrected  onc's  in  i'ig.  VII  .'(.a),  i.ource  d  i  spl  aceinont  spectra  are 
obl.t'ned  siinpl'"  b'.-  correcting  ali  the  wav  hack  to  t.he  source,  v.'ith  riit- 
t  o]  lotting  parameters  invoivi-d; 

-  -inpv  '  •  S  •  it  r  ‘  ^ 

wliti’i.  P  exp(.eK  vi.j  j  ane .t<^  t  ic  attenuation,  G  Ls  geometrical 
spreading  as  delinetl  above,  S  is  radiation  pattorn  coeificif'nt  (G.b) 
times  free-surtace  aiip.  I  i  t  i  ca  t  ion  t2.0)  divided  by  a  possible  VL.tori,, 
partitioning  of  f  nei'g;.  (J7).  v  is  v;ave  '.-elocitv,  />  is;  density,  and  U  . 
is  t  .le  ob.served  displacement  spectiaun.  In  applying,  these  co  t  n  <■  t  inns, 
'..'ith  a  .Mociel  1  t  t.f-ir  1.1 1  ion ,  we  g,et  source  d  i.sp  1  aci'inent  spectra  a.s  .slio'  .. 
in  Fig.  VI1.4.J,  lyhcre  a  soismic  moment  of  the  order  of  Mm  (  llh  '' 

dyne 'em)  is  indicated.  In  using  i  lu‘  Hanks  and  Kananiorl  (19/'*)  moment 
magnitude  re  1  a  t  i  onsli  i  [' 

•^w  “  •'  h  Mn  -  6.0 

wo  then  get  ).  ),  whii  li  is  quite  consistent  with  the  otiicr 

magqiitude  (stiniitcs  di.scussod  above.  For  .an  earl  liquake  this  s  i  r.e  ,  the 


Rl'.me  source  model  gives  s  corner  frequency  at  0.8  }\z  for  100-1'or 
stfess  drop. 

It  is  not  possible  from  Fig.  Vll.'t.S,  however,  to  deterinltie  corner 
frequency  with  any  reasonable  accuracy.  The  reason  for  this  is  partly 
low  frequency  noi.se  (as  mentioned  above),  but  prin:nrlly  the  fact  tha . 
the  Q-model  (Q  120  •  f^-^)  most  probably  is  not  applic.Thle  feu: 
frequencies  below  1  Hz.  A  lower  limit  in  Q,  possibly  even  combined  with 
an  increase  towards  lower  frequencies  (Aki,  1980),  v\;ould  yield  the  lot- 
frequency  asymptotic  effects  called  for  by  the  coinihonly  icceptcd  source 
models.  The  sensitivities  and  the  uncer tainties  involved  here  are 
properly  illustrated  by  the  low  frequency  differences  between  the  tto 
correction  models  in  Fig.  VII. A. 4:  one  order  ol  magnitude  difference  -r 
a  distance  of  about  300  km  (Molde,  Sulen)  and  two  orders  of  mugnitud- 
dii  ferences  at  about  twice  that  distance  (NORF.SS)  . 

Anotther  question  that  it  raised  from  the  present  observ.uions  is 
concerned  with  the  small  differeitces  between  the  hr  amplitudes  at 
Molde/Sulen  as  compared  to  NORKS.S .  In  the  corrected  spee'!-i  {me  Fly, 
VTl.4.5)  this  shows  up  in  the  higher  NORESS  iovehs.  in  spite  ot  ' he 
fact  that  the  same  time  window  has  been  used  in  the  two  cases  i as 
compared  to  using  comparable  group  velocity  window-;  .Sin-co  we  have  : 
found  any  techtiical  reasons  for  this  differeiu;e  (s-i-.di  as.  errors  in 
gain)  we  assume  that  the  rea.son  must  be  tied  to  Lg  ■  a-.-e  p'opaga "  i  on 
characteristics  that  are  not  being  adequately  predicted  i--  tiu-  ode', 
used  here. 


Concluding  remarks 

The  questions  raised  here  cill  for  cotitimu-d  (ttort.^;  aiiiuti  it.  icsol'.-h- 
c.xi.sting  uncert.-iintie.s  in  our  knowl  et'.g.e  about,  am- 1  .as  t  i  a  ;  t  >  nua  t  i  on 
over  n  wider  range  of  f requonc ie.s  . 
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Fig,,  VII  ■  4  ■  4  ■  Ground  motion  displacement  (in  m'sec)  vs.  frequency  foi 
the  same  data  a.s  shown  in  Fig.  VII. 4. 3  (but  for  the  Z  compovient  only^  . 
path  corrected  hack  to  a  reference  di.stance  of  10  km  from  the  source 
using  the  tollowing  two  0  models;  (1)  Kvainme  and  Havskov  (1988)  and  (, 
Sereno  et  al  (1**88).  The  straight  lines  correspond  to  a  slope  propor¬ 
tional  to  . 
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Fle. ■  VII  ■  4 . 5  ■  Source  displacement  spectra  for  the  Aui’ust  8.  1988, 
earthquake,  with  corrections  for  geometrical  spreading  (Herrmann  and 
Kijko,  1983)  and  anelastic  attenuation  (Kvamme  and  llav5;kov,  1988).  Tl.e 
plot  indicates  a  seismic  moment  of  the  order  of  10'^  Newton  -  me ters 
(equivalent  to  10^“^  dyne-cin)  .  The  straight  line  indicates  a  slope 
proportional  to  ,  and  curves  1  through  h  are  Moldc  ,  SuKn,  NORESS  it’ 
and  NORESS  HF.  respectively. 
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ARCESS/riNF  SA  .-ut  a  vs 

This  contribution  coniprisf's  t  v;o  n(!parat;c  i  nvcsl  i  j^aci  ons  related  Lcp 
analysis  of  events  ro.'orded  on  the  three  regional  arrays  NORESS .  AHeKSS 
and  FINESA  in  Fennoscandia .  The  first  investigation  is  an  evaluation  ni 
the  performance  of  tin-  recently  \ipgraded  FINESA  array  in  Fin!  arid, 
whereas  the  second  investigation  utilizes  data  recorded  s  imultaneousl'.' 
on  all  three  arrays  in.  producing  joint  event  locations. 


An  evaluation  of  the  perfortnance  of  the  upgraded  FINESA  array 

A  description  or  tlui  KINES.A  .irrav  is  given  in  Korhonen  et  nl  (1987),  '  i’ 

early  1988,  the  ,-coii'cury  of  the  FINESA  array  was  expanded  hy  atiding 
five  elements  to  !  1)0  .rrrav.  as  sliown  In  Fig.  Vll.b.l,  The  FINESA  ari'.y 
geometry  curronL'v  copprises  15  vertical  only  se  i  sniome  ters  within  <)n 
.■‘.’pc’-r  ure  of  2  kr  . 

FINESA  data  are  recorded  on  magnetic  tape  at  the  annjy  .site,  and  the 
rape  recoiding  is  normally  event  triggered  by  a  built-in  voting 
detfctor.  In  order,  liowevt.  r .  to  I'fcpei'lv  (-v.ali  .to  the  performance  of 
the  upgraded  FINESA  array,  data,  ■..cro  recordcil  cent  inunusly  for  a  la-d.v 
period  during  Mai'ch  .S-.M  .i|  I'h’.t;  1 1,^.  tapes  weii'  ]>l.ivod  hack  and 

checked  at  NORS.AR ,  and  approx  i  in.'i ;  c  Iv  a5%  ot  tin-  d.ata  for  tl.i.s  li-d,;V 
period  could  be  iccovered  and  Wtu.'  hence  '.nhjected  to  detection 
processing.  The  remaining  of  t.hc  dat  a  could  not  ho  read  due  to 

various  prt.i  lems  wifli  ‘ho  '..ipt.u;,  liku-  parity  errors,  etc. 

A  beam  deployment  comprising,  !'?  In  ari.s  f  OG  colieiu  iit  ,  G  incoherent)  was 
used  for  the  detect  ion  procf  ss  i  nr,  ol  the  continuous  FINESA  data.  The 
beam  deployment  used  i  .s  in  ag,ref’iiieMt  with  the  recommendations  hv  Kv.rn 
et  al  (  1987)  and  Kv.Trna  (I'tgS).  The  deteetiou  processing  results  in 
list.s  with  attributes  for  each  detected  signal,  like-  detection  time, 
sign.i]  frequoiicv  phase  veloeitv  and  arrival  azimuth,  rheso  lists  woi. 
compare-d  ag.ainsf  the  rcp,ional  Finnish  huLletin,  it.sned  by  the  llni\'er- 


sity  of  Helsinki,  and  the  results  of  the  comparison  are  given  in  Tablt 
VII  .  5.1.  Only  those  bulletin  events  occurring  when  th.e  FiNESA  array  v;a 
operating  properly  are  included  in  the  t.ible. 

Signals  detected  on  FINESA  were  a.ssoc  i ;  ed  to  the  Helsinki  bulletin 
evi  nts  by  requiring  a  reasonable  match  of  FINESA  d.  ,ect:nn  par.inieter.^ 
tarrival  time,  phase  typo  from  velocity,  .ind  arriv.il  .i:;  i  ii..n  h '/  ..  i  t  h  tii' 
■or  responding  ones  predicted  frtm  the  i  nf  oriiitii  i  on  in  t![t  Helsinki 
bulletin.  From  Table  VII. ‘'.I  we  see  that  out  of  thi‘  101  .'t  tt-rence 
evet^ts  listed,  99  had  at  least  oik-  d('te(t«.-d  1’-  or  S-phase,  i.e.,  96  jit- 
>.'ent.  T'JO  of  the  tour  events  th.it  were  iiot  di- ■  t  ec! .  occurred  at  t  he 
Lahnaslampi  mine  in  Finland  (6ti  2‘’N,  2f.0‘'’E,),  at  a  '.ii.stance  of  322  kir 
!  rom  FINESA.  Most  blasts  at  Lahints  1  am()  i  art>  ipiite  .srra  1  1  ,  and  are  not 
detected  by  FINESA.  The  t.wo  remtiining  e'.’eiit.s  ./<.  re  both  simll  ones 
'iiiagnitvide  les.s  th  in  2  iOr  one  i  veiit  ;  ii.igniiiidi  I'l.j-  riv.  :.  to;  •  n 

•  at  rangt  s  mor-  '  h.iii  'X'  l.m  t  ro!i  j-' I  VF.s.y . 

ilnse  resnl‘'.s  for  th.,  vtpr.t  .■vied  1  lN|-:;.\  .■<  .|in  •  .,•!■■  .ui.-.  '.  . 

i  :r  '  s;  ig.it  ion  based  on  lo  o]  i  n.ii.e-:  i-lN'lS.'-  In  i  :  i-'' 

■  h  'he  o  r  i  r,  i  n.'i  1  a, 'ray  gt.oinetrvi  i-iuic  i -.di  d  'li.i'.  a.  i-er  •  v'  ■■  ■!.• 

event.s  list.;!  in  the  Helsinki  iH.:ii<..tin  wi- .  .  tUt.-  t  L-<.i  k--.  t  le. 
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Joint-  event  locations  from  three-array  data 


Data  recorded  at  FINESA  were  used  together  with  NORKSS  and  ARCESS  dat  i 
in  cissessing  tlie  capahi  li  ri  es  of  this  three-array  networV;  in  locating 
events  in  the  Fe'.nno.scandian  region.  A  set  of  10  events,  for  which  there 
was  at  least  one  detected  phase  for  each  array,  was  selected  for  an 
event  location  experiment .  The  events  are  listed  in  Table  VII.  5. 2  and 
shown  in  Fig.  V 1  I  . 1 . 2 .  The  event  magnitudes  range  from  less  than  2.0  to 
3,2,  The  origin  times  and  geographical  coordinates  for  the  10  cents 
are  taken  from  rlie  If  lsinki  hnlliMin. 

The  continuous  r  i  ocei-.s  inr.  ol  data  recorded  at  o.'ich  of  the  three 
regional  -U'rav.s  in  Fenino.sc.nifii  a  provides  e.s  t  i  ina  l.c.s  of  arrival  times  c.id 
back  aticriths,  n- ir.tnie  i  e  rs  together  with  the  a.s.sociated  uncertain- 

:  ie,-  wc-re  used  a.'--  fhi  T'iVJll.Of:  program  eh-ve  1  opc-ti  t)y  Bratt  ,n;d 

Baci'.e  (I'^PS'i.  T'!',-.;;!  i  ncinporat  cs  the  arrival  time  and  aaimuth  d..t.'i 
into  a  g.euera  1  i::c(i  -  i  :r.  i-rse  Location  estimation  .scheme,  aiici  can  be 
c-.ppiiid  Ic  both  "  I  -  arr-iv  and  :rn  1 1  i  p  1  e  •  .-i  rra  v  data. 

Fig.-,.  VI:  b,3a,  Fo  and  Ic  siiow  AR'IESS.  NTtRF.RS  .atid  FIN'ESA  data, 
re  .s  g'ec  t  i  ve  1 V  .  :  ,.  :<i  3  i  r.  Tabic  Vll.'.'.2.  The  |i<nie  1  s  show  on  i.lie  top 

•itre  •  P-wavi  hra.ip.';  ioi  thr.  dii'i  rent  freepienc',  haitd.s  .  Tltc  beans  we  r. 
■ste.-red  a'- o  rel  i  T.g,  to  t  pc  pi.  ,s<-  it-,  .m-l  ..ritirath  nf  the-  peak  ot  t!,- 

1  rt  -  .'i].ie'C  r  ri.in  comp-at.t-d  ,:c;  p.if  f>(  tin-  oitliite-  detect  i  ein  proi-estniig  Th.- 
t  h  ft  ■  1  o'We  L'  t  t'.ic  e  -  <:  o  i  t  h'-  p.'O'i  ■  1  to;  ,  i  ,  i . ,  1  <  ,i  I  ' .  - 

lor  titre'e  diifefi-nt  tre-qtien.  v  h.n.els  .  Tiu-.  oote-c'  ion  time-s  fot  ti,-- 
pha-Sf-s  usod  in  •  liu-at  io-.  e-‘-.  pc  r  i -.se  i  ;  itr>  iiuifr-  ic,-  ,i  ro-,.,-  ,, 
fiptires  .show  that  thi.-.  -.-v--!  •  is  •  a  i  i.r-l.  d  ,.  i  •  b  ,  i.ig.b  Si!  •  •- 

elosc'St  eiiTei-,-  (  1  1  N  i-.:!.-',  .  '.-.•lu  I't  .c.  it  is  ::  ■  p.,  t -.wc,  -th  :  nto-. 

1 1.  i  .s  de  t  e-  c  t  e  d  elue*  to  the-  t, .  .'R  g .  i  i  i  -  tb.it  is  ,  n  ;  i  d  t  1 ; e>n  t- 1 ,  !  u  <  i  v  o  a:  ■ 
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Table  VII.  5. 2  gives  the  results  of  the  location  exp.  rimciit  .  On  ttie 
-ivcrage,  the  joint  three-array  locations  deviate  from  the  network 
locations  published  in  the  Helsinki  bulletin  by  16  km.  Two-array  and 
one-array  locations  were  computed  for  all  combinations  of  events  and 
array  sub  -  networks ,  also  using  the  TTAZLOC  algorithm.  Thi  resulting 
average  deviations  from  tlie  r.etwork  solutions  arc  26  and  68  km. 
respectively . 

The  results  for  onc-.arrav  and  two-arrav  locations  .'.re  Ln  general 
agreement  with  what  has  previously  been  reported  (e.g.,  Mvkkeltvei C 
and  Ringdal  (1988)  found  an  average  deviation  of  34  km  fiom  the 
riel.sinki  bulletin  loi-ations,  using  data  irom  seven  regional  events 
recorded  at  NORESS  and  AKCF.SSi.  The  improvement  in  the  locatior. 
accuracy  when  invoking  datr.  I'rom  three  arrays  is  significant,  and  we 
otisider  the  result.s  reported  here  as  quite  promising,  when  taking  ti.r 
tdllow'ing  into  account:  T!;.  ai'tival  timt's  used  wt-re  those  dc  t  r;uined 
lUt  omatically  by  tlie  onlin.'  pri'ce.s.s  i  ns . .  It  is  cone..  ; ’.'ah  i  e  th.a'..  liuman 
i  ;U' '■  r  vent  ion  for  adiustment  ol  ariiv.il  time.s  f’l.ment  of  th. 

j'jto/tatlc  procedure  would  iinpru’.’i  t  lu-  loc.ition  .-st  i  ,i’a; '-s .  Only  standard 
l  avol  time  tahles  for  the  pliases  Pn.  Sn  and  i.g  were  u.st  <1  I'he 
Mit  1  odiic  t  ion  of  region,!  I  i  .:<-d  travel  time  tables  is  'ii:<  to  I'lsult  ir. 
iir.prov‘‘ ments  .  i'innll’,.  master  I't  loe.ition  .<;chem(  s  of  ■.  iiiou-,  kind’ 
i.old  con.s  i  de  rab  le  proniis.  and  .ire  eypiTtefi  to  iurtluf  <-nh,in,'e  tiK 
e.ip’.iiii  1  j  t  i  e.s  of  accui'itelv  1  oe;i  t  i  ni;  regional  •'■..■nt:., 

S .  Mvrl •  It-  :  : 

.1.  Kver 

;  .  Kvmi  na 

M,  It:  i  .  fn  .  v  .  b.  !  :  t  ; 


PeffUjuicf’.'. 


Bratl  ,  8.P.  nut  '1' .  t:  .  B.aelu-  .l-iB8i;  I.oe.e  i  ng  e'.'i  nt  s  with  spaisi 

lU'twork  of  rer.ional  event:..  Bull  .  Seism.  Sue.  A.i'.  i  ,  S  ,  /'80-  "'t' , 

.''oiTionen,  H.,  S.  I'irhonen,  I-' ,  H  i  ngda  1  ,  S,  Mykkellveit.  T.  Kvm'rn.i  .  i' .  i-.' . 
h.irsen  and  R.  P.'cil.sen  i  1  Vfi  ’  )  :  The  •  .iie.sa  arrav  u'd  le  1  i  in  i  :’a  r  v 
r. 'suits  of  d.ifa  .-ui.'!  I  v.s  '  s  .  I'n  i  .  <>l  Hi-l.sinki.  1  ::t  .  e'  Se  i  si:'<Ti  og  ■  ■  , 
i'>  po r(  S  •  1  ■>  . 


Kvavna,  T.  (198S):  On  exploitation  of  sinall-aporCure  NORESS  type  arviv,- 
for  enhanced  P-wave  dftectabillty .  Semiannual  Tech.  Sumin.,  1  Ocr 
1987  -  31  Mar  1938,  MORSAR  Sci.  Rep.  2-8;/88,  Kjeller,  Nor-.'ny. 

Kvsina,  T.,  S.  in.bsgaard,  S.  MyUkeltveit  and  F.  Ringdal  (1987):  Toward 
an  optimum  be;im  di;p  1  oyment  for  NORESS;  experiment.s  with  a  North 
Sea/Uostcan  Norway  data  base.  Semiannual  Tech.  Summ .  ,  1  Apr  -  3(i 
Sep  1987,  NORSAR  Sci.  Rep.  1-87/88,  Kjeller,  Norway 

Mykkeltveit,  S.  and  F.  Rlngdal  (1988):  New  re.sults  trom  pi'oce.ss  i  np,  of 
data  recorded  at  t  b.c  new  ARCESS  regional  array.  Semiannual  Tecii. 
Summ.,  I  Oct  198:'  -  31  Mar  1988,  NORSAR  Sci.  Rep.  2-87/88, 
Kjeller,  Norw.ay. 


Date 

Time 

Lat. 

(ON) 

Lon . 
(°E) 

Magn. 

Dist.  P 

(km) 

-det.  S-det, 

38/03/08 

06.01.42 

61.2 

27.1 

<2 

61 

X 

88/03/08 

07.50.37 

62.8 

29.1 

<2 

215 

X  X 

88/03/08 

13.03.50 

62.2 

23.3 

<2 

169 

X  X 

88/03/08 

13.29.02 

64 . 2 

28.0 

<2 

322 

- 

88/03/08 

14.25.10 

62.1 

26.4 

<2 

7  5 

X  M 

88/03/08 

15.02.39 

60.3 

24 . 8 

<2 

1  4  5 

M  X 

88/03/09 

12.05.43 

64 . 2 

28.0 

<2 

yA2 

88/03/10 

08.47.04 

59 . 3 

2  7.2 

<2 

247 

X 

88/03/10 

09.22.08 

59 . 3 

2  ; .  6 

<2 

253 

!*’ 

88/03/10 

09.57.40 

59 , 2 

27.6 

<2 

264 

X  o’. 

88/03/10 

10.41.14 

59 . 3 

27.6 

253 

y. 

88/03/10 

11.07. 37 

59.2 

2  7.6 

2  1 

264 

X 

88/03/10 

11.14.23 

59.3 

27.6 

<2 

253 

>;  X 

38/03/10 

11.25.34 

59.5 

25.0 

2 

224 

X  >; 

88/03/10 

11.49.54 

59.3 

27.6 

253 

X  X 

88/03/10 

12.05.20 

59  .  5 

26.5 

■)  •> 

218 

88/03/10 

12.07.05 

61.2 

2,8.9 

- 

154 

X  X 

88/03/10 

12 . 10.50 

59.3 

28  .  1 

2.  1 

264 

X  X 

88/03/10 

16.03.30 

64 . 3 

24.0 

<2 

335 

X 

88/03/10 

16.20.56 

62.0 

24 . 4 

.  ;> 

108 

X  M 

88/03/10 

18.16,  15 

65.8 

.  9 

.  '•) 

4  9] 

88/03/10 

18.29.30 

6 .  1 

20 . 6 

2 

88/03/10 

20.27.28 

63.6 

26 . 2 

<2 

2  '.ii 

88/03/11 

08.18.59 

62 . 9 

2  5 . 9 

: ..  5 

X  X 

88/03/11 

09.23.17 

67.6 

54 .8 

■ 

s  * 

88/03/11 

09.24,26 

67,6 

5.  0 

- 

•83 

88/03/11 

09.25.40 

67.6 

34 . 0 

2  .  , 

783 

X  X 

88/03/11 

09.48,06 

61,4 

34 . 3 

2  .  5 

-.39 

X  x 

83/03/11 

10.21,09 

62.2 

25.9 

•22 

8  5 

X  X 

88/03/11 

10.21.35 

59.3 

2  7.6 

2  2 

.>5  3 

X 

88/03/11 

10.56,54 

59  ,  5 

2  5 . 0 

2  .  1 

2  24 

X  X 

Sa/03/11 

11.27,25 

59.3 

2  7.6 

<2 

.-5  5 

V 

88/03/11 

11.46 . 58 

69.4 

30,8 

2  .  3 

9  13 

X 

88/03/11 

12.03, 37 

6  3.2 

2  7.8 

2.  5 

.'  1 

88/03/1 1 

12.17.09 

59  .  5 

25.0 

1 

1 

88/03/11 

12.33.24 

60 . 8 

29  .  3 

•  >  ) 

!  .5  8 

88/03/1 1 

12.57. 59 

59  .  3 

27  6 

•"2 

.'5  3 

•8/03/11 

13.3  5.05 

59  ,  5 

28  .  1 

2 

26  . 
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ion  pi 
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1  in:  x.-\  tLa..  I  <i: 

the  period  8-21  Mar 

1988.  fhe  t.ihle 

1  i  .s  (  .s 

tlie  10  5  tv. 

.11  f  s  ( » (  [  he 

He  1 s i nk i 

bullet  ill  tha'  occurred  while  llie 

KINKS  A  .s  v  . 

t  cm  v:<t.s  ope I'a  1  i i 

properly 

during  the 

14 -day  period. 

llie  distance  1  roir 

the  FINESA  arix 

is  piven 

for  each  event.  The 

table 

iridicate.s  whetht  i' 

of  T>or  a  F  or 

phasf'  v;as 

de  t  ec  ted 

on  FINESA 

,  that 
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as5.,o<'  i  u «  -  i 

x'itli  *lu’  e\cMt 

fj'  u';.  t  ion. 
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Date 


Timo  L.at  .  Lon.  Mctp.n.  Dist.  P-riot..  S-dot. 

(°N)  (^F.)  (km) 


88/03/12 

09.12. 28 

59.4 

28/4 

2.1 

261 

X 

X 

88/03/12 

09 . 59 . 59 

6  <4  .  7 

30 . 7 

2.9 

4  51 

X 

X 

88/03/12 

10.43.17 

..9 . 5 

25,0 

2.2 

224 

X 

X 

88/03/12 

10.48 . 21 

61.8 

36.1 

2.4 

533 

X 

X 

88/03/12 

11.03,58 

!)8  .  1 

3  3.2 

2  ,  1 

8  i  5 

X 

X 

88/03/12 

11.11.12 

n8  .  1 

4  4  4 

4  4  . 

2 . 9 

8  1  5 

X 

X 

88/03/12 

12 .25.01 

-0 . 3 

4  7  3 

2 . 3 

24  7 

X 

X 

88/03/12 

12.40.45 

(;  4 . 6 

'30,6 

718 

- 

- 

88/03/12 

12.41.07 

‘  \b 

30 . 5 

.  4 

16 

X 

X 

88/03/12 

14.  ill  38 

■'  7  .  1 

20 . 6 

<.2 

684 

X 

38/03/13 

06.49  10 

hi  .1 

3  3.7 

2.4 

7  86 

x 

88/03/14 

09.01.  .1 

jO  ,  3 

2  7. 6 

v  7 

253 

X 

X 

88/03/14 

09 . 13  21 

62 . 8 

22.!. 

<2 

2  36 

X 

X 

88/0  3 /14V: 

09 . 22  .  16 

50  ,  3 

2  7.2 

2  .  1 

24  7 

X 

X 

88/03/14 

10.32.41 

>9 . 3 

2  7.0 

<2 

253 

X 

X 

88/03/14 

10.35  25 

3 

i  7 . 6 

<2 

2  53 

X 

X 

88/03/14 

12.  41,  . O' 

vO,  6 

.30.0 

2.1 

298 

X 

X 

38/03/14 

13  .  10  ,  5,1 

59.5 

2  5 . 0 

2.6 

224 

V 

X 

88/03/14 

14.07.27 

/*.3 

28.  L 

<2 

264 

V 

v; 

88/03/15 

08.59.57 

.  6 

34 . 0 

<2 

783 

- 

X 

88/03/15 

09.06. 42 

6  7.6 

34.1) 

2 . 4 

783 

- 

X 

88/03/15 

10.31.30 

:;9 . 2 

7 . 4 

2  .  5 

260 

X 

83. 0  1/1  5 

11.34  36 

v9. 5 

2  o  .  4 

2.  5 

218 

X 

X 

88/03/15 

11.41 . 57 

bO .  5 

2  5 . 0 

•> 

106 

X 

X 

88/03/15 

11.26. 30 

01.6 

21  .  , 

.l  34 

X 

X 

88/03/15 

12.10,  38 

■'0,2 

2  7  .  5 

• .  / 

204 

X 

86,  03/15 

12.11.  :  7 

>0 . 2 

,1  7  ^ 

,,  > 

26-4 

X 

88/03/15 

12.  19.15 

,  3 

,'7  7  ,  , 

•5  > 

7 

X 

88/03/15 

12.33,  55 

.  5 

2  1  .  7 

;  "‘■H 

X 

N 

88/03/1 5 

12,40.  '30 

50.7, 

28  .  '■ 

> 

V 

88/03/ 15 

13.  16  .  s5 

■  1  ,  3 

. 

X 

88/03/15 

13.22  .  14 

n  .  9 

3(1 , 6 

?  .  0 

;*  4‘t 

\ 

X 

88,'03/15 

13.51.30 

5  5  ,  o 

'.>5  K 

•> 

; 

X 

X 

88/03/15 

14.0'' .  3  5 

'0 . 5 

5  5  .  •  1 

)  > , 

X 

88/03/15 

14,20. 59 

60. 0 

4  0  2 

'i 

,  i 

X 

88/03/15 

14,36, 30 

6  3  .  1 

"'2 

/ 

>  j 

\ 

88/03/15 

14 . 39  .  3  5 

5  0 , 5 

7  i 

> 

•  1 

.  /  » 

S' 

88/03/15 

17.57.51 

6  5,8 

7  4  ,  , 

- 

\- 

88/03/16 

08.37,  1  5 

S') ,  3 

/  /  ■'  ■ 

) 

2*^ . 

X 

\ 

88/03/16 

719 .44 . 50 

f.O  .  (. 

j  .  '4 

X 

- 

*  The  orij; in  lime  foi 

’.his 

i  .’eni  i 

;  ir  i '  p  r 

\  [\l  <iS 

on  'Mi  1  f, 

in  rl 

Helsinki  b\il It- Lin 
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Table  VII  ■ D  .  1  . 


(I’.ir,e  .  ol 


«) 


Date 

Time 

Lat . 

Lon . 

Magn  . 

Dist . 

P  -del. 

S-det 

(ON) 

("E) 

(km) 

88/03/16 

10.25.03 

59.3 

27.6 

<2 

253 

X 

88/03/16 

10.45.^0 

60.9 

26.8 

<2 

72 

X 

X 

88/03/16 

11.26.47 

59.2 

2  7.6 

<2 

264 

V 

X 

88/03/16 

11.45.36 

63.2 

27.8 

2.5 

215 

X 

X 

88/03/16 

11.49.51 

59.5 

25.0 

<2 

224 

X 

V 

88/03/16 

23.04.26 

67 . 8 

20.0 

<2 

765 

- 

- 

88/03/17 

09.07.13 

58.3 

10.9 

2.7 

91/ 

X 

y 

88/03/17 

10.21.17 

69.6 

29.9 

2.0 

926 

X 

X 

88/03/17 

10.27.20 

59 . 2 

2  7 . 6 

2.3 

264 

X 

V 

88/03/17 

10.46.21 

59 . 2 

2  7.6 

<2 

264 

X 

V 

88/03/17 

11.18.48 

59.3 

27.2 

2.3 

.’4  7 

X 

V 

88/03/17 

12.02.23 

64.2 

28.0 

<2 

322 

- 

88/03/17 

12.02. 36 

59.4 

28.5 

2.1 

.■’64 

X 

88/03/17 

18.58.07 

59  ,  ■’ 

5.6 

3.2 

1135 

X 

X 

88/03/18 

05.16.20 

6  0  .  ..> 

34 . 7 

2.6 

9  52 

X 

X 

88/03/19 

10.04.08 

61.1 

30.2 

<2 

224 

X 

38/03/19 

10.05.02 

59 . 3 

2  7.2 

- 

24  7 

A 

X 

88/03/19 

12,15.34 

68.1 

33,2 

■ 

815 

X 

88/03/19 

12.15,39 

68.1 

33.2 

815 

X 

X 

88/03/19 

12 .39.09 

68  .  1 

33.2 

<2 

815 

X 

X 

88/03/19 

13.03.39 

67.6 

30,5 

- 

718 

N' 

\\ 

88/03/19 

13,03.54 

6  7.6 

30 . 5 

- 

718 

:-;8/03/19 

13.07.00 

61.5 

30.6 

2.6 

.3  a  5 

X 

88/03/19 

13  42.33 

67 . 6 

30 . 5 

<2 

718 

- 

88/03/20 

04. .*5.  17 

67.7 

33.  / 

2.5 

786 
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Fig..  VI  I.  5,?.  Th( 
NORESS,  ARCF.SS  inc 
in  thf  TTAZLOr,  1  o( 


Event  3  at  ARCESS 


Be  oi 


Bqqi 


AOZ 


1988-075/1^ :2?:50.000 
ARCESS 


BU  BP  2 .0-  1  4  .□  NKPL  5  PAS  il 


BU  BP  5.0-  5.0  NPOL  5  PAS  I 


BU  BP  5.0-  10.0  NPOL  5  PAS  I 


SB  0.  -  t.S  MPOL  2  PAS  I 


BU  BP  1.5-  ,5.5, 


.  u  5 .0-,  1,0 .0  NPOL  5  PAS  ,1 


SECO^CS 


KIr,  VlI.S.3a.  ARCESS  data  foi'  event.  1  in  T.ihli  V I  [  . .  '  The  p.inel 
shows  on  top  three  P-beam.s  steered  towards  the  epieenter.  Lor  tliree 
di  fferent  filter  bands.  Tlie  three  bottom  tra.  »  s  correspond  to  tliree 
different  filters  applied  to  the  vertical  sen.sor  at  site  .AO.  The 
tletection  times  (by  the  automatic  online  proce.ssor)  ioi'  the  phases  Pi 
and  Lp,  are  indicated  by  arrow.s. 


Event  3  at  NORESS 


B  €E?  Q  I 


B  e  o  t 


B  Q  o  m 


eu  BP  S.O-  5.0  NPOL  i  PAS  I 


id 


BU  BP  5.0-  10.0  NPOL  5  PAS  ( 


13  10 
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Event  3  at  FINESA 
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i  M'sA  d  , :  ;  .  The  s  L ;  a' 
■\l  '  t  Ki  .  VI  I  Ml  .  n 


VI I  ■  6 _ Co  111  pa  ra  1. 1  vi  anal  vs  i  s of  N()!.SAR  and  ( ■  la  1  eiiha  r  Lt'  ma  p.n  i  rude  ^ 

for  Siiar.an  Kiver _ o  xijJ  os  j  oiis 

Introduction 

The  seismic  Lp,  cav-  p :  tc's  in  the  I'ont  i  iienta  1  1  i  t  hosplie  La-  and  eai. 

Pe  observi'i  a';  f  if  iv.'av  as  lOlK'l  |■.l^  in  sliiel'l  and  stLiliK'  piatfai'm  are.is 
(NuLtli,  1^73:  B.uuLa'.a.  rdt  ,  1983).  I.;’,  i  .s  p,enii;aliy  eons  i  dt-rt-d  'o  Consist 

of  a  supe' rpos  i  t  i  on  of  many  h  i  f.he  r  -  mocie  surfacr  v/nvt  s  oi  pj-oiip  velo¬ 
cities  near  3.0  kin  -  .s  ,  .and  i  t  .s  rad'.-uiun  i.s  th>-refort  o\]'ecl.ed  to  be 
more  isotropii-  tli.m  tlsat  of  P  w.-n.-.-s.  Thu:.,  full  at’,  imntiial  (  (ivc-rape  is 
not  essen'i.ai  for  reliable  de  term  i  nat  i  on  of  hr,  m.-ipnitud.  .  Pnrtherirore  . 
Lf,  is  nor  atfee'  <1  liv  lati''.,!  he :  ,>  ro,r-,  ne  i  :  t  es  in  :h-.-  nppe;  mantle, 
•..•'li.n  can  prod'..','-  i  eig  i  oe  n.ss  i  ipg/de  f  oens.s  i  nj,  .  *!■  e  ;  s  oi  i' -  ■.■.'a'.a  • ,  .,!'d 
t  1k'  •  f o ft  c on :  r  i  h' : '  .  ■  a  :  :  e ;  ■ '  :  C ■  ,n:  mu'e  r '  .a  ■  ■  '  -  li.i  d  mp 

r : •  r '.\\:\  t"  h-  S  , 

;  r  n"  i  ( ii'  i  s  r .ih  1  (■  e';  ■  :  ,n ,  t  <  . > :  i-  ; . 

hi’',  a  . hi  Tt'  ■■  t  S  i  '  I  I  ,  1  i  ,  .i;  ■  P,  ■ 

S'  1'  i  pa  1  a  t  i  risk  <■  s.n  1  e s  .  ‘a i:  :  1  i 

r  o  .  -I  •  im. , ;  i  m,  ' .  ■  , i  •.  i . ■  1  d'l  of 

P  i  : a',1... 1 1  V  1  a 3  i  ’  I  i  I'  -m  -  t  : 


t  :  ;  fee  i  . ,  i  on  i !  :  a  .  .  -  i  m  : 

mt  !  i  lOfl  j :  1  ini'-,,  ' ,  '  i  j  .  .  ,  :  .  ;n  I "  ■ 

(  1  98  '  ,1  :  aid  Pi  nrda  1  ,  d  h  '■  n  i  '  ' 
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have  assessed  the  ett'ects  of  introducing,  station  coriections  for 
individual  array  elements  and  cpicentral  distance  corrections  in  the 
estimation  procedure.  The  precision  in  the  e.stimates  has  been  invest  l- 
■’,ated  takinp,  into  account  the  s  i^nal  -  to- noise  ratios,  and  a  comparattv^e 
analysis  of  NORSAR  and  Grafenberg  Lg  measurements  h.is  been  carried  out. 


Data  sources 

The  NORSAR  array  (Bunguni,  Husebye  and  Ringdal,  1971)  was  established  in 
1970,  and  originally  comprised  22  subarrays,  dtployed  over  an  area  of 
100  km  diameter.  Since  1976  the  number  of  operational  subarray.s  has 
been  7,  comprising  altog, ether  92  vc-rtical -component  SP  sensors  (type 
HS-10).  In  this  paper,  analysis  has  been  restricteci  to  data  from  these 
7  subarravs.  Sampling  rate  for  the  NORSAR  SP  data  is  20  samples  per 
second,  and  ell  data  are  recorded  on  digital  magnetic  tape. 

The  Grafenberg  array  (Harjes  and  Seidl ,  1978)  was  established  in  1976. 
and  today  comprises  13  broadband  seismometer  site.s.  tbri.  e  of  which  a''e 
3-component  sy.stcms.  The  instrument  response,  is  f  l.r.  to  velooitv  froit 
uiiout  20  second  period  to  5  Hz.  Sampling  rate  is  20  .samples  per  second 
,'.nd  the  dat.a  are  recorded  on  digital  magn.eiic  taptv 

The  location  of  NORSAR  and  .,.nberg  relative  to  S>  m  i  p.i  la  t  i  nsk  is 

shown  in  Pig.  VII.  6.1,  whett.:  al.so  tiie  propagation  p..ths  to  the-  two 
.irrays  are  indicated. 

Ba.se-d  on  1  SG  and  NEIC  reports,  a  toteil  of  99  >  vi  nt  s  .  pre).S’.iincd  t  ■'  be- 
nuclear  explosions  at  the  Shagan  River  area,  have  been  .selected  as  a 
liata  base.  The-  time  .span  is  from  1969  to  Se-pteunber  19,  l'»88,  when  the 

second  Joint  Verification  Experiment  (JVE)  explosicn  w.t.s  carried  out. 
T.ab  1  e  Vll.6.1  lists  the  d.ites  of  therse  events  toge-the-r  with  pertinent 
measuremeuits  discussed  leitei  In  the  text. 


Data  analysis 

All  ."Vdilable  x\‘coi'div.gs  from  NORSAK  atid  CRF  have  bt'>n  ar.;ilyzcd  for  the 
event  set  of  94  Shagar  River  explosions,  using  tlie  procedure  described 
by  Ringdal  and  Hoklnnd  (1987). 

Briefly,  this  p'ocedr.re  compri.sts  filtering  all  array  cliannels  witli  a 
0 . 6  -  3 . 0  H?,  bandpass  filter,  computing  R.MS  value  of  each  filtered  trace 
in  a  2-minute  fg  window  (starting  1?  inin  after  P  onset  for  NORSAR,  14 
min  for  GRF)  .  and  coirnens.it  i  ng  for  background  noise  preceding  P-onset. 
The  1-g  magnitude  is  i  nen  e.stiniated  bv  logarithmic  averaging  acros.s  each 
1  rrc  ■.  . 

The  total  uumije  i  oi  a.'.’ailable  recordings  witli  .sufficient  signal  -to- 
noise  ratio  to  allow  reliable  Lg  measurement  war:  70  for  NORSAR 
(starting  in  l,9.’l)  atul  AO  f(,)r  CRF  (starting  in  1976). 

vsTiilo  the  NORSAR  avniv  cont  i  gura  i  i  on  has  been  stable  over  the  time 
period  considered,  rl^.  C;PF  array  initiallv  conpristd  only  the  four 
instruments  Al.  -  A-l  .  .n<i  was  latei  expanded  to  iti;  full  conf  i  gnra  t  i  on 
of  13  sites.  In  ordc!'  to  reduce  as  far  as  pos''.ihle  the  bias  due  to 
changing  array  t.iud  i  gurnt  i  ons ,  we  have  therefore  computed  station 
coriectioiis  for  e/ten  iudividual  (’Rf  ^uue'or  (Table*  VM  1.6.2)  and  tipplicd 
tho.se  in  ehe  arr.iv  .av.  raging,  procedure.  A  similar  .';ot  of  correction.s 
for  NORSAR  are  i.i.stei.i  in  'i'.mle  VII. 6. 3.  In  practice,  the  introduction 
of  station  cor  rt  c  t  ion.*,  h.’.s  made  little  difference  ft):'  the  NORSAR 
magnitude  estimate.s,  hut  had  a  .significant  effect  for  CRF. 

The  eifects  of  t'pici.’nt  .  a  1  il  i  .s  (  tnuu*  >'  !  fer>  lu't-.s  on  t  lu*  I.g  ir.a  r.n  i  tude 
estimates  have  al.so  hi  'n  ;i.*:.se*>';ei! .  The  diat.inie  I'ur  rf-e  t  iuci  B(A';  is 
determined  through  iNe’iJi.  ]986h): 

■’.h'.)  f  .s  i  re  •,/ 1  ’  1  )  /  i  ni  A()/l  1  1  )  I ■  exp  [  y  ( i’.q  )  ] 

Af)  is  ;  iu*  '.iistanec  (I’m  to  a  fixed  reference  loe.it  ion  within  tlu- 
epicentra'  ,irea  (foi  Srm  i  pa  1  a  t  i  nsk  we  have  used  SO‘’N  ,  49"K)  and  is 
the  d  i  *  ■  ,,:  .  I*  ( Inr  .  <  ,, 


:  iie  event  . 


y  is  the  coefficient  of  .anelasti*' 


attenuation.  We  have  used  7  =  0.001  km'^,  which  is  near  the  value 
obtained  by  Nuttli  (1986b)  for  1  second  Lg  waves  for  paths  from 
Semipalatinsk  to  Scandinavian  stations.  Note  that  a  very  accurate  value 
of  7  is  not  required  when  considering  a  limited  source  region,  as  the 
efiects  of  small  variations  in  this  parameter  on  tiie  res'.ilting  mg(Lg) 
values  are  negligible. 

The  Lg  magnitudes  at  KORSAR  and  CRF  of  events  in  the  data  base  are 
listed  in  Table  VII. 6.1.  Since  thc^se  estimates  take  into  account  both 
station  terms  and  epicentral  distance  corrections,  they  are  slightly 
different  from  values  published  earl  if)-,  but  nevertheles.s  in  good 
agreement . 

Table  VII. 6.1  also  contains  estimated  standard  deviations  of  the  Lg 
magnitudes,  taking  into  account  both  Che  scattering  across  each  arrav, 
tlie  signal  -  to-noi  se  ratios  anci  the  variance  reduction  obtained  by  the 
averaging  procedure  'see  Appendix).  We  empha';ir.e  that  tluse  standard 
devLationr,  are  indicati\'e  only  of  the  precision  c>t  measurement,  and 
should  not  be  interpreted  as  being  reprosental ivc  oi  the  accuracy  oi 
these  magnitudes  as  .source  size  estimators.  We  notf  that  magni  t  vides  oi 
the  larger  explosions  may  be  measured  with  very  high  precision,  v.-h.ertas 
the  uncertainty  is  greater  for  the  smaller  events,  due  tn  the  lower 
signal- to-noise  ratios.  It  is  also  cle.ir  that  the  NtlRSAR -based 
estimates  are  mere  precise  than  those  using  GRF  data,  especially  for 
events  for  which  full  GRF  tirray  recordings  are  not  .ivti  i  l.ihl  c . 

Fig.  VIT.6.2  shows  a  scatter  plot  of  NORSAR  versu.s  '..RF  iiei c,n i  ti.'.de.s  toi 
all  common  event.s.  The  strt.ight  line  represents  a  l<,ist  squares  fit  to 
the  data,  assuming  no  errors  in  NORSAR  magnitvide.s .  'a\  note  that  the 
two  arrays  show  excellent  const  .st.ency ,  although  thort'  i  .s  .iorae  increase 
in  the  scattering  at  low  magnitudes.  The  standard  diviatiou  of  the 
differences  relative  to  the  least  squares  fit  is  O.OA5  magnitude 
ttnit.s.  Also  there  is  no  significant  separation  between  events  from  NF 
and  SL'  Shagan  with  regard  to  the  relative  Lg  magnitudes  observed  at  the 
two  arrays. 


In  Fig.  VII. 6. 3  a  similar  plot  is  shown,  including  only  "wel 1  -  recorded" 
events,  i.e.,  requiring  at  least  5  operational  GRF  channels  and  a 
standard  deviat  ion  ot  each  array  estimate  not  exceeding  0.04  niagnitucie 
units.  The  slope  of  the  .straight  line  fit  has  been  restricted  to  the 
same  value  (1.13)  a.s  in  Fi,g.  VII. 6. 2.  We  note  that  there  is  a  sig¬ 
nificant  reduction  in  the  sc.ittei-,  and  the  standard  deviation  of  the 
residuals  is  only  0.032  magnitude  units.  Thus  the  i.g  magnitudes 
measured  at  the  two  .tt rays  show  excellent  consistency  for  high  signal - 
to-noise  ratio  events. 

The  slope  (1,15)  oi  'lie  straight-line  fit  in  Fig,.  VII, 6. 2  is  slightly 
gre.-iter  thrin  1 .  (iO ,  a  tender^cy  also  noted  by  Ringdal  and  Fyen  (1988): 

The  interpretation  of  this  observation  is  somewhat  uncertain;  a 
possible  explanation  is  scaling  differences  in  the  Lg  source  spectruir 
(Kv*rna  and  Ringdal,  1988).  in  combination  with  the  response  dif¬ 
ferences  of  the  'T'KSaR  and  GRF  instruments.  We  have  attempted  to 
compare  the  two  data  sets  -aftor  .".ijusting  the  GRF  recording.s  to  a 
\'0RS4R-type  response.  However,  the  results  were  inconclusive  since  th.’ 
GRF  s  igna  I  -  to  -  no  use  ratio  iheii  became  too  low  for  the  smaller  events. 

Fig.  VII. 6. 4  illustrates  the  pattern  of  P-Lg  bias  in  the  Shagan  R  i '.’e  > 
area,  using  m^,  values  computed  at  Blarknost  fMarshalL,  jiersonal 
communication)  together  with  combined  NORSAR/ilRK  Lg  magnitudes.  The 
latter  have  been  dei  ived  by  ad  justing  the  GRF  magnitude  s;  to  an 
"equivalent"  NORSAR  v.tlue  using  tlie  .st  ra  ight  - 1  i  ne  relation  of  Fig. 

VII.  6.  3,  and  then  calculating  a  weighted  average  using  the  inverse 
variances  (Table  Vll.A.l)  as  weighting  factors.  Fig,  Yll.6.4  inc.lvules 
all  events  of  >  5.6,  assuming,  ^ithl■r  two-array  olisevvat  i  ons  or 

very  precise  Lg  iiK'a.sori  ineiits  t  roin  one  array  (<'  0.0'''). 

Although  both  the  mg  ''allies  and  the  l.r.  m.agni tude.s  have  been  revised 
relative  to  those  used  in  earlier  studies,  lig.  Vtl.6.4  confirms  the 
observations  previous;  ly  made  regarding  the  systematic  difference 
between  P-Lg  re.sidual.s  from  NE  and  SW  Shagan.  In  the  NF.  area,  mgfPt  ’ " 
goncrall'y  lower  ’  ii.u'  ni(Lg,),  whereas  t  lie  opposite  behavior  i  .s  .seen  in 
the  SW  portion.  The  .IVE  explosion  of  14  September  1488  has  a  P-Lg  bias 


of  0.06  which  is  clo.se  to  the  aver.-ige  for  tlie  SW  region.  Fur therinore , 
there  appears  to  be  a  transition  zone  between  the  two  portions  of  th.: 
resc  site,  where  tlie  residuals  are  close  to  r.ero. 

Conclus ions 

From  this  and  previous  studies,  we  can  conclude  that  the  bg  RMS 
estimation  methods  provide  very  stable,  mutually  consistent  results 
vdien  applied  to  two  widely  separated  arrays  (NORSAR  and  (dlF)  .  Tlii  s  i.s 
of  clear  significance  regarding  tiie  potential  use  of  such  Lg  measure¬ 
ments  for  yield  estimation.  Further  research  will  be  directed  toward 
expiinding  the  data  base  by  conducting  similar  studie.s  using  other 
available  station  data  as  well  as  studying  Lg  recordings  from  ocher 
test  sites.  In  particular,  seismic  data  that  might  become  available 
from  USSR  statioirs  in  the  future  would  be  of  importance  both  in  furtiici 
assessing  the  stability  of  the  estimates  and  to  obtain  l.g,  magni  tvides 
for  explosions  of  Ic'v;  yields. 

F.  Ringdal 

J .  Fven 
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I'abl  e  VI I  ■  6  ■  1  ■  List  oi  presumed  explosions  at  tlie  SIkii’.ii  River  test 
■  Mca  neeir  Sem  Lpalat  insk ,  USSR.  The  m|j  values  are  tliose,  published  in  the 
iSC  buJLetins  for  events  prior  lo  1986,  and  are  otlierwise  taken  trom 
NEIC/PDE  reports.  NORSAR  and  Grafenberg  Lg  RMS  magnitudi-.s  are  given  (or 
all  events  with  available  recordings  of  sufficient  s i gna L - to  -  no i se 
ratio.  The  number  of  data  channels  used  and  tlu-  estimated  precision  c.f 
measurements  (see  Appendix)  are  given  for  eacli  magnitude  value.  (Pagv  1 
of  :-)  . 


No . 

ORIGIN 

DATE 

ORIGIN 

TIME 

MB 

****  NOHSAR  **** 
M(LG)  N  STD 

■k  it  -k  If  it 

M(LG) 

GRF 

N 

***** 

STD 

51 

04/22/81 

1 

17 

11 

6.0 

5.907 

28 

0.022 

5.956 

11 

0.027 

52 

05/27/81 

3 

58 

12 

5.5 

5.456 

27 

0.023 

- 

- 

- 

53 

09/13/81 

2 

17 

18 

6.1 

6.114 

29 

0.008 

6 . 109 

9 

0.015 

54 

10/18/81 

3 

57 

2 

6.  1 

5.984 

34 

0.010 

5.956 

9 

0.021 

55 

11/29/81 

3 

35 

8 

5.7 

5.545 

28 

0.121 

5.512 

12 

0.192 

56 

12/27/81 

3 

43 

14 

6.2 

6.071 

34 

0.009 

6.050 

9 

0.021 

57 

04/25/82 

3 

23 

5 

6 . 1 

6.078 

35 

0.008 

6.069 

10 

0.017 

58 

07/04/82 

1 

17 

14 

6.1 

- 

- 

- 

- 

- 

- 

59 

08/31/82 

1 

31 

0 

5.3 

- 

- 

- 

- 

- 

- 

60 

12/05/82 

3 

37 

12 

6  .  1 

5.988 

31 

0.019 

6.001 

13 

0.020 

61 

12/26/82 

3 

35 

14 

5 . 7 

5.655 

39 

0.080 

5.598 

13 

0.067 

62 

06/12/83 

2 

36 

43 

6 . 1 

6.073 

25 

0.009 

- 

- 

- 

63 

10/06/83 

1 

47 

6 

6.0 

5.867 

19 

0 . 033 

5.851 

11 

0.040 

64 

10/26/83 

1 

55 

4 

6  .  1 

5.999 

33 

0.021 

6.035 

11 

0.020 

65 

11/20/83 

3 

27 

4 

5.5 

- 

- 

- 

- 

- 

- 

66 

02/19/84 

3 

57 

3 

5.9 

5.723 

29 

0.038 

- 

- 

- 

67 

03/07/84 

2 

3  9 

6 

5.7 

5.695 

29 

0.065 

5.575 

12 

0.108 

68 

03/29/84 

5 

19 

8 

5.9 

5.899 

29 

0.012 

5.961 

13 

0.043 

69 

04/25/84 

1 

9 

3 

6.0 

5.869 

35 

0.008 

5.804 

13 

0.031 

70 

05/26/84 

3 

13 

12 

6.0 

6.073 

33 

0.007 

6.132 

13 

0.015 

71 

07/14/84 

1 

9 

10 

6.2 

6.055 

32 

0.007 

6.066 

12 

0.015 

72 

09/15/84 

6 

15 

10 

4  .  7 

- 

- 

- 

- 

_ 

- 

73 

10/27/84 

1 

50 

10 

6.2 

6.082 

33 

0.011 

6.143 

13 

0.016 

74 

12/02/84 

3 

19 

6 

5.8 

5.881 

29 

0.020 

5.864 

12 

0.036 

75 

12/16/84 

3 

55 

2 

6.1 

6.046 

29 

0.010 

6.037 

13 

0.014 

76 

12/28/84 

3 

50 

10 

6.0 

5.982 

35 

0.009 

5.944 

13 

0.021 

77 

02/10/85 

3 

27 

7 

5.9 

5.801 

40 

0.024 

5.800 

13 

0.058 

78 

04/25/85 

0 

57 

6 

5.9 

5.859 

29 

0.045 

5.848 

7 

0.047 

79 

06/15/85 

0 

57 

0 

6.0 

5.976 

30 

0.009 

6.031 

13 

0.017 

80 

06/30/85 

2 

39 

2 

6.0 

5.928 

30 

C.009 

5.905 

12 

0.016 

81 

07/20/85 

0 

53 

14 

5.9 

5.858 

37 

0.013 

5 . 867 

12 

0.031 

82 

03/12/87 

1 

57 

17 

5.5 

5.215 

33 

(-.076 

- 

- 

- 

83 

04/03/87 

1 

1  7 

8 

6.2 

6.051 

33 

0.008 

6  .  126 

11 

0.017 

84 

04/17/87 

1 

3 

4 

6.0 

5.898 

3  3 

0.020 

5.912 

12 

0.026 

85 

06/20/87 

0 

53 

4 

6  .  1 

5.968 

36 

0.007 

5.943 

10 

0.028 

86 

08/02/87 

0 

58 

6 

5.9 

- 

- 

- 

5.856 

11 

0.022 

87 

11/15/87 

3 

31 

6 

6.0 

5.973 

37 

0.008 

5.983 

13 

0.022 

88 

12/13/87 

3 

21 

4 

6.1 

6.091 

31 

0.010 

6.066 

12 

0.015 

89 

12/27/87 

3 

5 

4 

6.1 

6.046 

31 

0.011 

6.032 

13 

0.019 

90 

02/13/88 

3 

5 

5 

6 . 1 

6.042 

26 

0.009 

6.047 

13 

0.029 

91 

04/03/88 

1 

33 

5 

6.  1 

6.067 

31 

0.007 

6.076 

13 

0.014 

92 

05/04/88 

0 

57 

6 

6.1 

6.040 

31 

0.008 

6.064 

13 

0.020 

93 

06/14/88 

2 

27 

6 

4.9 

- 

- 

- 

- 

- 

- 

94 

09/14/88 

4 

0 

0 

G  .  0 

5.969 

37 

0 . 0 1  n 

5 .970 

12 

0,04' 

Table  VI I  ,  ]  .  rp,i_('_(.  ,>  o(  p  ; 


ANNEL 

NO 

BIAS 

N 

STD 

1 

0.15 

24 

0.029 

2 

0.15 

31 

0.031 

3 

0.19 

24 

0.042 

4 

0,08 

19 

0.034 

5 

0.01 

12 

0.046 

6 

-0 . 11 

18 

0.030 

7 

0,01 

16 

0.041 

8 

0.09 

15 

0.036 

9 

-0.09 

19 

0.039 

10 

-0.15 

13 

0.024 

11 

-0.04 

7 

0.033 

12 

-0.17 

12 

0.039 

13 

-0,12 

14 

0.045 

Tab  Lci  VI I  ■  6  ■  2  .  List  of  .station  tcrin.s  (.station  KM.S  I.,;  va  Uio  minus  ar;.nv 
average)  for  the  (Jrafenberg  array.  The  13  individual  vertical  component: 
■sei  .smometers  are  Listed  in  the  sequence  Al-^.  Ul-'j  ntd  C 1  -  ^  .  Tlio  bia.s 
values  are  ba.sed  on  high  signa  I  - 1  o- no  i  se  ratio  event  .s  recorded  by  at 
Ica.st  10  cb.annels.  The  number  of  observation;:  and  the  .sample  standard 
deviation  i.s  listed  for  each  instrument. 


CHANNEL 

NO 

BIAS 

N 

STD 

1 

0. 

05 

49 

0.051 

2 

0. 

11 

31 

0.044 

3 

0. 

17 

23 

0.044 

4 

0. 

04 

6 

0.025 

5 

0. 

10 

45 

0.028 

6 

-0. 

01 

49 

0.060 

7 

0. 

01 

43 

0.033 

8 

0  . 

08 

39 

0.045 

9 

-0. 

01 

42 

0.029 

10 

0  . 

13 

42 

0.036 

11 

0. 

00 

34 

0.047 

12 

0. 

05 

43 

0.042 

13 

0. 

03 

48 

0.038 

14 

-0. 

11 

48 

0.028 

15 

-0. 

01 

49 

0.033 

16 

-0. 

01 

49 

0.035 

17 

0. 

03 

49 

0.032 

18 

-0. 

02 

49 

0.040 

19 

-0. 

02 

43 

0.033 

20 

-0. 

01 

44 

0.043 

21 

-0. 

05 

45 

0.034 

22 

-0. 

05 

45 

0.024 

23 

-0. 

03 

44 

0.049 

24 

-0. 

04 

46 

0.022 

25 

-0. 

10 

45 

0.031 

26 

0. 

02 

45 

0.037 

27 

-0  . 

07 

44 

0.027 

28 

-0. 

08 

45 

0.023 

29 

-0  . 

02 

45 

0.031 

30 

-0  . 

02 

45 

0.038 

31 

-0. 

06 

42 

0.031 

32 

-0. 

Oa 

42 

0.025 

33 

-0. 

03 

41 

0.047 

34 

-0  . 

02 

43 

0.033 

35 

-0. 

04 

44 

0.029 

36 

0. 

01 

40 

0.055 

37 

-0. 

04 

21 

0.031 

38 

-0. 

05 

32 

0.030 

39 

0. 

01 

20 

0.064 

40 

-0  . 

01 

19 

0.046 

41 

0  . 

02 

18 

0.036 

42 

0  . 

05 

20 

0.029 

Tahl<;  VI  I.  6. 3.  i.i.st  <,j  sl,ili,.n  tfiiiis  (lit  a  lion  RMS  wiKu  iTiiniia 
average;)  for  I  lie  N'ORSAR  ari'ay,  Tlii-  a?  individual  so  israoire  t  v  r  <;  ai  o 
listed  in  'he  staiideird  .sequenc't;  I  r  I'avs  0]A  t  hroujvli  06C)  .  i'he  ii  i 

viilues  are  based  on  ei’ents  with  high  signal  -  to- noi.se  r.itio  ( i,g 
magnitude  >  5.8;.  The  number  of  observ.it  i  ons  .'uid  the  sample  siandai 
deviation  are  li';:ed  for  each  i  i\s  (  ruiiu  ii  f  . 


EAST  LONGITUDE  (DEGI 


V T  I  .  -j  .  J  .  I-or;U  i  on  of  Llii'  NORSAK  .ind  Ur;)  t  i  iibi  ;  .  ■■:)' 
))o  St  inip.i  1 . 1 1  in.'.k  test  sit  e. 


i;i  rflc'it 


MAGNITUDE  COMPARISON 

SHACAN  RIVER  EVENTS 
5=  1  15  I-— 0.901  Z>=  0.045  S3 


Fig.  VII  .  h  .  2  ■  riot  nl  i  :)-,i  f  (CKF)  vt.T.stis  .NdRSAK  (NAO)  l.g 

magnitudes  tor  Sliagan  River  exp  I  i  on.s .  Tin  figure  inriutUs  ,ill  eoir.inoi, 
events  in  Table  VI I  .  b  .  i  ,  Event  .s  in  the  NF.  and  SW  parts  of  Shagan  a  i\ 
marked  as  filled  squares  and  opet.  squares,  respeo  t  i  ve  1  v .  Tlie  sLraiglit 
line  (slope  1. .  1  b )  repr<-sents  a  liast  .squarcss  fit  to  the  data,  assumiu; 
no  error  in  NORSAR  Lg  measurement ss .  The  standard  deviat  tor.  ■'!  t  lu 
residuals  along  the  vert  leal  axis  relative  to  t  hi'  straight  line  is 
O.Oab,  and  the  dotted  lities  <'orre‘.ponfi  to  p  I  us /in  i  nus  two  s;andar<l 


GRF  M(LG) 


ISC  MB 

SHAGAN  RIVER  EVENTS 
BIAS  RELATIVE  TO  MAO/CRF  LG 


Fl£ ,  VI  1.6.4.  Ploi  ol  P-Lg  m.'igni  Muio  s  (ISC.  maximum  1  iko  1  i  1  uu.d 

minus  NORSAR/Grafeiibrrg  Lg  magnilmk>s)  as  a  iunrlioii  of  e\'f'nt  local  ion 
(Marshall,  personal  communication)  within  the  Shagan  River  area. 

Plusses  and  circles  correspond  to  residuals  greater  or  less  than  the' 
average,  respectively,  with  symbol  size  proportional  to  the  deviation 
All  events  of  m^(Lg)  >  5.6  for  which  we  have  precise  locations  liavc 
I  been  included,  assuming  eitiier  two-array  observations  or  very  precise 

Lg  measurements  from  one  array.  Tlie  JVE  explosion  is  especially  marked. 
Note  the  systematic  variation  from  NE  to  SW  Slia;;an.  v/ith  an  apparent 
transition  zone  in  between. 


Aopendlx  to  Section  \/II.6 


Ip  this  appendix  we  develop  an  approxiiiuit.e  expression  i'oi'  the  uncer¬ 
tainty  in  the  RMS  Lg  magnitude  estimates  described  earlier.  We  first 
consider  the  case  of  a  single  sensor  mea5-.urement ,  and  afterwards 
address  the  array  averaging  procetture. 

Denote  by  X]^(tl  the  recorded  signal  in  the  "Lg  window",  and  assume  that 
this  is  composed  of  .a  noise  component  x^(C)  and  a  signal  component 
x;jit)  as  follows: 

xi(t)  ■-  x^Ct)  +  X3(c)  (1) 

Hero,  we  assume  that  the  noise  component  x^It)  can  bt'  monel  led  ,i.s  a 
tero-mean  random  process  whicli  is  stationary  over  a  Lime  interval  long 
enough  to  include  both  the  Lg  window  and  a  suitable  noise  window 
preceding  the  P  >nset .  The  signal  X3(t)  is  considered  a  zero-mean 
random  process  defined  in  the  Lg  time  window,  atid  t>c-ing  itncorrel  ated 
wi  th  X2(t’)  . 

We  can  thus  obtain  an  estimate  of  the  mean  square  '’alue  .\3  of  \3(t)  h\- 

X  j  -  X-)  -  Xp  (?  ) 

where  X3  is  the  mean  square  value  of  xj^tt)  in  t  lie  .sigr,al  winder.,  and  X) 
is  the  mean  square  value  o[  xylt)  in  the  noise  window. 

The  Lg  RMS  magnitude  is  then  (apart  from  an  tidditi'-  const, mt  ) 
determined  as  log^Q 

We  now  make  the  assumption  that  the  quantities  ( 1-- I . ’)  each 

follow  a  lognormal  distribution,  when  consid<-rc-d  as  random  vari.ibles. 

We  emphasize  that  this  assumption,  which  is  reasonat)  1 e  in  view  of 
empirical  studies  of  logarithmic,  amplitude  patterns  of  signals  and 
noise,  represents  an  approximation  only.  Thus,  we  know  th.)t  the 


difference  botv^en  t\vo  lof'noi-mal  vnriaV)Les  is  visnany  not  ai\other 
lognormal  variable.,  but  for  our  purposes  this  approximation  is  useftil  . 

We  tray  thu';  wriin  (using  natural  logarithms); 

logXj^  is  i  ^  1 . 3  (3) 

Note  that  using  'ta >0  the  variance  of  logXj^  corresponds  to  o 
representing  the  vari.ince  of  the.  log  RMS  estimate. 

The  mean  and  Verin.nccs  ol  (  !ie  ro.spectivt-  vatiahl  e.s  can  then  be 
expressed  by  (Aitchi''on  and  Brf)wn.  1969); 


EXj^  -  e  t  1 

i  “• 

J  ,  .  .  ,3 

(^) 

0 

var  Xi_  -  (hXj^)'^  ■ 

.? 

'  -  1)  i  - 

1 . 3 

(5) 

F  rom 

eq.  (?)  We  fur  t  lie 

rmort- 

obta i n 

EX3  =  F.X|  F.X; 

(6) 

var  X,  • 

VcJf  i 

io 

t?) 

Coirb 

ining  .'  a)  a;id  ( /' )  . 

this 

Kad.s  to  the 

r(  Lat  i on ; 

(EX2  -  E.X’i)-  ■  -  1)  -  (EXi)-'  ■  (e'"’''  -  I)  + 

(KX^)-'  •  -n  (8) 


Substituting  EX^  and  KX2  by  i  ho  i.bserveti  values  X|  and  Xp ,  respec¬ 
tively,  and  assuming  .small  values  of  Oj  (i  ■  1 . 3)  we  obtain  from 

(8)  the  following  simplified  relation; 


105 


^3 


2 


A  ^  O 

(Xi  -  X2)2 


(9) 


which  represents  an  approximate  expression  for  the  variance  of  logTx^ 
Note  that  (9)  is  developed  using  natural  logarithms,  it  applies  without 
change  if  base  10  logaritlims  are  used  throughout. 

Although  we  have  used  a  number  of  simpl i f icat.ions  in  arriving  at  (9;, 
simulation  experiments  using  randomly  generated  distributions  have 
shown  that  this  formula  gives  a  useful  approximation  to  the  actual 
scatter  in  the  estimates  within  a  reasonable  range  of  parameter  values. 

^  A. 

We  note  that  in  case.s  of  high  signal-to-noise  ratio.s,  (i.e.,  X]^  »  X)), 
we  obtain  from  (9)  a-^  »  thus  the  noise  variance  has  no  sig¬ 

nificant  effect  on  the  Lg  magnitude  variance.  On  t!ie  other  hand,  as  the 
signal-to-noiso  ratio  becomes  small,  the  variance  ej2  will  increase 
rapidly . 


In  rhe  array  averaging  procedure,  we  assume  that  tlu'  leii;:  Wj2  jg 
reduced  in  proportion  to  the  number  of  array  element  s.  wIi.Toas  we 
consider  to  represent  mainly  a  systematic  noist  fluctuation  that  is 

not  reduced  through  array  averaging. 

A  A 

Defining  the  signal-to-noise  ratio  a  by  a  ■-  Xj/X;,  ,.i;d  di.ioting  1,'y  N 
the  number  of  array  elements,  we  thus  olttaLn  i  ivim  >  '  '< 


(<72^2  ■  +  £^22 


{n  -  !)■ 


(10) 


As  a  numerical  example,  consider  the  JVE  explosion  (.event  99  in  Table 
VII. 6.1)  . 

For  NORSAP ,  we  have  etitimated  a  -  13.12,  with  N  --  37,  and  we  assume 

0.04,  02  -  0.08.  Formula  (10)  thttn  gives  03  -  o.OlO. 


i06 


For  GRF,  we  have  q  -  3.03,  with  N  12,  and  the  same  input  a  values  as 
above  then  give  “  0,043.  Thus,  the  estimated  uncertainty  of  the-  CPF 
Lg  magnitude  is  considerably  greater  than  that  of  NORSAR,  the  main 
reason  being  the  lower  signal-to-noise  ratio  for  GRF. 

Reference 

Aitchison,  J.  and  J.A.O,  Brown  (1969):  The  Lognormal  Distribution, 
Cambridge  I'ni  vt  rsity  Press,  UK. 


(if  1  hf  strniglit  line  lias  been  restricted  to  the  value  (ibt, lined  in  IM 
V  1  I  .  f) .  .  Note  that  tlie  scatter  in  the  data  ha;  beer,  r;  i  gii  i  f  i  rniil  I  v 
red, iced,  and  the  r.tandard  deviation  in  the  vertical  direct  inn  is  on  I 
(l.0!2  magnitude  units  for  this  data  set. 


average,  respocLively ,  with  symbol  size  proportional  to  the  deviatim; 
All  events  of  mp,(,Lg)  >  5.6  for  which  we  liave  precise  locations  lo-nr 
been  included,  assuming  either  two-array  ob.servat  Ions  or  very  precise 
Lg  measurements  from  one  array.  Tlic  JVE  explosion  is  espec  ially  marked. 
Note  the  systematic  variation  from  NE  to  SU  Sha.gan,  -..'Ltli  an  apparc-nt 
transition  zone  in  between. 


follow  a  logiioirraal  distribution,  when  considi-rod  as  random  vari.iljle 
We  emphasize  that  this  assumption,  which  Is  ri-.ison;.()Ie  in  view  of 
empirical  studies  of  logarithmic  amplitude  patterns  of  signals  and 
noise,  represents  an  approximation  only.  Thus,  we  know  that  the 


For  NORSAR ,  we  have  estimated  q  -  13.12.  with  N  -  37,  and  we  assume 
^  0.04,  (72  =“  0.08.  Formula  (10)  then  gives  03  ==  ^hOlO. 


